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A  modern  computer  based  automatic  data  acquisition/con¬ 
trol  system  was  installed  at  the  Department  of  Aeronautics* 
Combustion  Laboratory.  This  system  utilizes  an  HP-85  desktop 
computer  as  system  controller  for  the  HP-3054A  data  acquisi¬ 
tion  system.  These  major  control  applications  include: 

(1)  particle  size  determination  by  measurement  of  scattered 
laser  light,  (2)  a  vitiated  air  heater  system,  and  (3)  a 
solid  fuel  ramjet  test  unit.  Details  involving  specific  data 
acquisition/control  techniques  for  these  applications  are 
discussed. 
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I .  INTRODUCTION 


This  thesis  project  was  undertaken  to  develop  and  install 
a  modern,  computer  based  experiment  control/data  acquisition 
system  at  the  Naval  Postgraduate  School  Aeronautics  Depart¬ 
ment,  Combustion  Laboratory.  Previously,  all  laboratory 
experiments  were  manually  controlled  and  all  test  data  were 
reduced  manually  by  hand  with  no  real  time  data  processing 
possible.  Precise  timing  of  test  equipment  and  control 
features  were  difficult  or  impossible. 

The  introduction  of  the  HP- 3 054 A  Automatic  Data  Acquisi¬ 
tion  system  and  the  integration  of  all  associated  test  equip¬ 
ment  to  a  desktop  computer  have  revolutionized  current 
operations  at  the  Combustion  Laboratory. 

Large  quantities  of  high  quality  data  can  be  acquired  at 
high  data  rates  and  stored  for  future  use  or  processed  real 
time  in  closed  loop  feedback  control  circuits  for  precise 
process  control  and  display. 

Other  electromechanical  equipment  such  as  Visicorders 
provide  analog  real  time  data  and  run  in  parallel  with  digi¬ 
tal  equipment  but  can  be  remotely  controlled  for  precise 
timing  of  experiments. 

Interactive  programs  on  the  HP-85  computer  allow  the  ex¬ 
perimenter  to  do  preliminary  calculations,  initial  experiment 
set-up,  and  to  perform  accurate  process  control  of  system 


functions.  On  line  data  reduction  and  display  of  results  is 
performed  by  the  system  CRT  display,  the  7225B  Plotter,  or 
the  system  printer. 

Four  basic  experiments/experimental  apparatuses  are  cur¬ 
rently  automated  using  the  automatic  data  acquisition  system; 
a  light  scattering  experiment,  a  vitiated  air  heater,  a  solid 
fuel  ramjet,  and  a  turbojet  combustor  test  rig.  The  first 
three  applications  are  discussed  in  subsequent  sections  along 
with  the  measurement  techniques  used  to  implement  the  process. 


II .  SYSTEM  AND  COMPONENT  DESCRIPTIONS 


To  produce  accurate  test  results  involving  high  speed 
machinery  and  processes  and  simultaneously  provide  precise 
timing  and  control  signals,  a  modern  data  acquisition  system 
has  been  installed  at  the  Combustion  Laboratory  at  the  Naval 
Postgraduate  School. 

Figure  1  shows  a  system  block  diagram  of  the  Hewlett- 
Packard  3054A  Digital  Data  Acquisition  System.  The  heart  of 
this  system  is  the  HP-85  microcomputer  functioning  as  system 
controller,  although  any  computer  compatible  with  the  Hewlett- 
Packard  Interface  Bus  (HPIB)  can  be  substituted  or  connected 
to  the  system  bus. 

A  description  of  each  equipment  connected  to  the  bus  is 
described  below. 

A.  3054A  DATA  ACQUISITION  AND  CONTROL  SYSTEM 

The  3054A  is  a  computer-based  automatic  data  acquisition 
and  control  system.  The  system  is  complete  with  a  system 
mainframe,  the  HP-3497A,  and  two  voltmeters  (the  HP-3456A 
and  HP-3437A) ,  that  are  interfaced  via  the  82937A  HP-IB  I/O 
card  to  the  powerful  HP  desk  top  computer,  the  85A. 

The  3054A  System  documentation  has  been  designed  to  co¬ 
ordinate  efficient  use  of  each  of  the  system  instruments, 
and  it  includes  a  specialized  software  support  package  that 


3054A  System  Block  Diagram. 


Figure  1 .  Data  Acquisition  System 
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is  designed  to  use  the  System  instruments  to  achieve  specific 
functions  and  operations. 


B.  HP-3456A  DIGITAL  VOLTMETER 

The  3456A  Digital  Voltmeter  is  a  3%  to  6%  digit,  inte¬ 
grating  voltmeter  with  relatively  high  speed,  10  parts  per 
million  (PPM)  basic  accuracy  and  100  nanovolt  sensitivity. 

Its  DC  and  True  RMS  AC  voltage  maximum  of  1000  volts  and  re¬ 
sistance  measurement  capability  reduce  the  amount  of  signal 
conditioning  necessary.  It  can  detect  100  nanovolt  changes 
in  100  millivolt  signals  at  speeds  of  48  readings  per  second. 

This  capability  is  required  for  measuring  thermocouples  with 
the  3054A  system  to  better  than  0.1°C  resolution.  By  se¬ 
lecting  4*j  digits  of  resolution,  reading  rate  is  increased  to 
330  readings  per  second.  The  3456A  has  selectable  integra¬ 
tion  times  from  0.01  to  100  power  line  cycles  and  provides 
input  guarding  to  give  140  db  of  common  mode  rejection,  pro¬ 
viding  for  accuracy  in  the  presence  of  noise.  Two  other  j 

i 

noise  reduction  features  are  provided  by  digital  averaging 
and  an  analog  input  filter. 

Up  to  350  readings  can  be  stored  in  the  3456A  Read/Write 
memory  in  addition  to  providing  a  programmed  time  delay  be¬ 
tween  readings.  A  "voltmeter  complete"  signal  provides  the 
necessary  means  to  synchronize  the  operation  of  the  3456A 
with  the  3497A  to  allow  analog  scanning  rates  up  to  300 
channels  per  second. 


C.  3437A  DIGITAL  VOLTMETER 


The  3437A  is  a  high  speed,  3^  digit,  DC  voltmeter  which 
provides  precisely  timed  sample  and  hold  readings.  It  can 
perform  a  variety  of  tasks  from  scanning  DC  inputs  to  AC 
waveform  analysis.  Repetitive  signals  with  frequency  com¬ 
ponents  up  to  1  MHz  and  low  frequency  transients  longer  than 
500  ysec  can  be  rapidly  digitized  and  analyzed.  Parameters 
such  as  RMS  value,  residual  DC,  harmonic  content,  and  peak 
values  can  be  readily  and  automatically  determined. 

Using  the  HP-85  controller  and  a  "Fast  Handshake"  trans¬ 
fer  mode  of  operation,  reading  rates  of  up  to  4000  readings 
per  second  are  possible  in  the  packed  binary  mode  of  operation. 

D.  3497A  DATA  ACQUISITION/CONTROL  UNIT 

The  3497  data  acquisition/control  unit  combines  precision 
measurement  capability  with  control  functions.  It  can  provide 
precision  measurements  of  strain  gauge  outputs,  thermocouples, 
pressure  transducers,  and  other  sensors,  and  can  also  provide 
digital  interfacing  and  control.  A  digital  clock  is  built 
into  the  3497A  mainframe  for  real  time  control  applications. 

Five  slots  are  available  for  plug-in  I/O  card  options  on 
the  3497A  mainframe.  Additionally,  the  3497A  capabilities 
may  be  expanded  through  the  addition  of  one  or  more  3498A  Ex¬ 
tenders.  Up  to  1000  analog  channels  and  1300  digital  channels 
can  be  obtained  with  the  addition  of  14  3498A  Extenders.  Each 
extender  has  10  slots  for  plug-in  I/O  card  options. 
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] 

The  3497A  mainframe  also  provides  several  key  features  ] 

that  make  it  very  suitable  for  use  in  computer  controlled 
data  acquisition  systems.  The  following  TTL  compatible  input 
and  output  terminals  are  standard  on  the  3497A:  external 
trigger  input,  external  increment  input,  break  before  make 
(BBM)  sync  pulse  (input /output) ,  voltmeter  complete  output 
signal,  channel  closed  output  signal,  and  a  timer  signal. 

Additionally,  the  3497A  has  multiple  programming  commands 
providing  versatility  in  process  control  applications. 

E.  349 8 A  EXTENDER  UNIT 

The  capacity  of  3497A  mainframe  may  be  expanded  through 
the  addition  of  from  one  to  fourteen  3498A  extender  units, 
each  having  10  slots  in  its  frame  for  I/O  modules  of  various 
types  to  meet  specific  application  requirements. 

Special  attention  to  the  3498A  extender  operating  manual 
is  required  when  adding  additional  I/O  cards.  Slot  numbers 
5  through  9  do  not  exist  for  digital  cards. 

F.  HP-8290 1M  FLEXIBLE  DISC  DRIVE 

The  mass  storage  requirements  for  the  system  are  handled 
by  the  8290 1M  dual-drive  disc  system  which  supports  two  5% 
inch  flexible  double- sided,  double  density  discs,  providing 
a  total  of  approximately  550K  bytes.  Drive  #0  is  established 
as  the  default  mass  storage  medium.  Over  220K  bytes  of  mass 
storage  are  also  available  in  the  HP-85  internal  casette  tape 
drive . 
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G.  HP-722 5B  PLOTTER 

The  HP-7225B  plotter  contains  a  17601A  personality  module 
and  is  therefore  a  microprocessor  based  HP-IB  plotter  that 
produces  high  quality  graphic  plots  on  any  size  chart  up  to 
210  x  297  mm. 

Pen  movement  is  as  small  as  .032  mm  for  high  resolution 
plotting.  The  plotter  instruction  set  includes  39  different 
instructions  to  provide  such  capabilities  as  point  digitizing, 
labeling,  character  sizing,  scaling,  and  window  plotting.  The 
7225  is  interfaced  through  the  HP-IB  providing  extreme 
versatility. 

Pen  velocity  is  programmable  from  10  mm/ sec  to  250  mm/ sec. 
Seven  different  dashed-line  formats,  symbol  mode  plotting  and 
user  defined  characters  aid  in  trace  identification. 

H.  HP-2631B  LINE  PRINTER 

The  2631B  is  a  180  character  per  second  printer  that  pro¬ 
duces  high  resolution  dot  matrix  printing.  Characters  are 
printed  from  one  of  two  128  ASCII  character  sets.  A  16  chan¬ 
nel  fixed  or  programmable  vertical  forms  control  (VFC)  is 
available  for  any  desired  page  and  text  length  with  several 
print  pitches  to  select  from.  An  underline  mode  is  also 
available. 

The  2631B  uses  a  bidirectional  printing  head  and,  addi¬ 
tionally,  monitors  incoming  data  to  determine  optimum  print 
direction  for  the  next  line.  Other  features  include  blank 


space  detection  and  suppression,  eight  print  densities  with 
underlining,  automatic  vertical  and  horizontal  tab  settings 
and  margins.  It  also  includes  an  end  of  line  (EOL)  wrap¬ 
around  feature.  When  the  printer  receives  data  which  would 
cause  printing  to  exceed  the  right  margin,  it  automatically 
moves  the  data  (which  would  ordinarily  be  lost)  to  the  next 
line. 

I.  HP-85A  COMPUTER 

The  heart  of  the  Data  Acquisition  system  is  the  85A  desk¬ 
top  computer  which  functions  as  system  controller.  The  com¬ 
puter  is  interfaced  with  system  components  via  the  82937 
HP- IB  I/O  card.  The  85A  uses  enhanced  BASIC  programming 
language  and  supports  the  PLOTTER/PRINTER  ROM,  INPUT/OUTPUT 
ROM,  MATRIX  ROM,  and  the  MASS  STORAGE  ROM,  providing  an  ex¬ 
tremely  versatile  instruction  set.  The  85A  also  features  a 
32x16  CRT  display  with  full  graphics  capability,  an  internal 
32  column  thermal  printer,  live  keyboard,  buffered  I/O,  pri¬ 
ority  intercept,  high  speed  tape  cartridge,  autostart,  and 
error  trapping. 

It  is  also  a  versitile  calculator  and  can  be  used  with 
virtually  no  previous  HP-85  experience.  The  user  can  halt  an 
executing  program  using  the  pause  key,  perform  a  series  of 
calculations  on  the  keyboard,  print  out  the  answer,  and  then 
resume  program  execution  by  pressing  the  CONTINUE  key. 

Reference  1  describes  the  HP-85  system  characteristics  and 
the  extensive  instruction  set. 


J.  HP-6942A  MULTIPROGRAMMER 

For  high  speed  data  acquisition  applications  the  HP-30 54A 
is  augmented  with  the  6942A  multiprogrammer.  The  6942A  is 
capable  of  functioning  independently  from  the  3054A  as  a 
master  control  unit  for  bidirectional  communication,  data 
processing,  and  data  transfer  between  an  HPIB  controller  (the 
HP-85  computer)  and  a  multiprogrammer  system.  It  can  be  used 
in  a  single-unit  system  employing  from  one  to  sixteen  plug-in 
I/O  cards,  or  in  a  multi-unit  system  consisting  of  one  6942A 
master  unit  and  up  to  seven  6943A  extender  units.  Each  ex¬ 
tender  unit  can  accommodate  up  to  sixteen  input /output  (I/O) 
cards,  allowing  a  multiprogrammer  system  to  be  expanded  to 
128  I/O  cards. 

The  6942A  employs  32  different  instructions  in  its  in¬ 
struction  set,  permitting  a  wide  variety  of  I/O  card  func¬ 
tions  under  program  control  of  an  HPIB  controller.  The 
instruction  set  includes  system  control,  output,  input,  card 
control,  and  system  timing  functions.  System  control  in¬ 
structions  establish  the  basic  operating  modes  of  the  multi¬ 
programmer  system.  Output  instructions  are  used  to  send  data 
to  output  cards  while  input  instructions  are  used  to  obtain 
data  from  input  cards.  The  card  control  instructions  allow 
access  to  various  circuits  on  the  I/O  cards  to  enable  close 
control  and  monitoring  of  their  various  operations.  The  sys¬ 
tem  timing  instructions  permit  the  sequencing  of  events  and 
the  measuring  of  elapsed  time.  A  detailed  description  of  the 
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instruction  set  is  described  in  the  6942A  User's  Guide  [Ref. 
2]  as  well  as  programming  details  for  the  HP-85  controller. 

In  addition  to  the  16  I/O  cards,  an  HP- IB  interface  board 
and  a  transmission  system  board  are  installed  in  the  rear  of 
the  6942A  unit.  The  6942A's  internal  mainframe  circuitry  in¬ 
cludes  a  CPU/ROM  board  and  a  control/RAM  board.  All  power 
supplies  for  the  I/O  card  slots  and  mainframe  circuitry  are 
built  into  the  6942A  unit.  The  CPU/ROM  board  includes  a  16 
bit  microprocessor  chip  which  decodes  and  executes  all  in¬ 
structions,  and  controls  all  data  transfers.  The  12K  x 16  bit 
ROM  (Read  Only  Memory)  contains  the  microprocessor  controlled 
programs  (firmware)  required  to  process  all  instructions. 

The  ROM  also  contains  diagnostic  programs  for  built-in  self 
test  and  signature  analysis/troubleshooting  capabilities. 

The  control/RAM  board  contains  a  2Kxl6  bit  RAM  (Random  Ac¬ 
cess  Memory)  which  is  used  for  the  temporary  storage  of  in¬ 
structions  and  data. 

Backplane  control  circuits  on  the  control/RAM  board  de¬ 
code  I/O  card  addresses  and  control  I/O  card  functions.  The 
HP-IB  interface  board  provides  the  bidirectional  communica¬ 
tions  interface  between  the  6942A  and  an  HP-IB  controller. 

Three  I/O  card  types  are  presently  installed  in  the  6942A 
multiprogrammer:  a  69736A  Timer/Pacer  card,  two  69751A  high 
speed  analog  to  digital  converter  cards  ,  and  two  69790B  4K  x 
16  bit  Random  Access  Memory  cards. 
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1.  69736A  Timer /Pacer  Card 

At  power-on  tha  timer/pacer  "wake a  up"  in  the  one 
shot  mode  which  produces  a  single  output  pulse  with  a 
crystal-controlled  duration  programmable  from  one  micro¬ 
second  to  about  eighteen  hours  (65535  seconds) .  Another 
mode  (the  continuous  mode)  can  also  be  programmed  in  which 
the  output  goes  high  for  the  programmed  time  interval,  then 
low  for  the  same  interval,  and  continuously  repeats  this 
cycle  until  the  card  is  reprogrammed  or  disabled  externally. 
The  card's  single  output  pulse  can  be  used  as  a  programmable 
delay  or  as  an  enable  signal  for  a  frequency  counter  card 
(69775A)  in  frequency  measurement  applications.  The  square 
wave  output  generated  in  the  continuous  mode  can  be  used  to 
pace  A/D  or  D/A  converters  as  they  analyze  or  generate  wave¬ 
forms  or  to  drive  frequency-to-voltage  converters  or  other 
devices  which  require  a  programmable- frequency  square  wave 
input . 

Figure  2  [Ref.  2]  shows  a  detailed  block  diagram  of 
the  timer  pacer  card.  In  the  Light  Scattering  experiment 
the  timer  pacer  card  is  programmed  to  the  continuous  mode 
and  drives  the  external  trigger  input  of  2  high  speed  A/D 
converter  cards.  The  timer  pacer  is  itself  externally  trig¬ 
gered  by  the  Reticon  Photodiode  array  #1  which  also  syn¬ 
chronizes  the  operation  of  Photodiode  array  #2. 
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Figure  2.  Multiprogrammer — Timer/Pacer 
[Ref.  2] 


2.  Analog  to  Digital  Converter  Card  I69751A 

The  69751A  is  a  12  bit,  analog  to  digital  (A/D)  con¬ 
verter  used  to  measure  bipolar  DC  voltages  in  one  of  four 
ranges:  ±100  MV,  ±1  V,  ±10  V,  or  ±100  V.  A  voltage  measure¬ 

ment  can  be  initiated  by  a  programmed  instruction  or  by  an 
external  trigger  signal.  Three  manual  range  switches  on  the 
card  select  either  the  ±100  MV,  ±1  V,  or  the  ±10  V  range. 

Input  voltages  in  the  ±100  V  range  are  connected  to  the  card's 
divide-by-ten  attenuator.  Sixteen  optional  DC  input  ranges 
(8  bipolar  and  8  unipolar)  are  also  available.  Instructions 
for  modifying  and  recalibrating  the  A/D  card  for  the  optional 
DC  input  ranges  are  described  in  section  11  of  Ref.  2. 

The  A/D  card  uses  a  high  performance  sample  and  hold 
amplifier  and  a  successive  approximation  A/D  converter  which 
provide  a  high  conversion  speed  and  excellent  accuracy  charac¬ 
teristics.  The  DC  input  voltage  is  guarded  to  maintain  high 
input  impedance.  Optically  coupled  isolators  are  used  to 
isolate  the  DC  input  voltage  from  the  data  lines. 

The  12  bit  digital  word,  indicating  the  magnitude  and 
sign  of  the  measured  voltage,  is  stored  on  the  card  and  is 
always  available  for  readback  to  the  multiprogrammer  main¬ 
frame  memory  using  card  subaddress  0.  As  many  as  1339  read¬ 
ings  can  be  stored  at  one  time  in  mainframe  memory.  The 
multiprogrammer  firmware  converts  the  digital  word  decimal 
form  for  readback  to  the  controller.  The  controller  can 


read  back  the  voltage  value  using  appropriate  HP-IB  extended 
talk  address. 

The  use  of  Extended  Talk  Addresses  and  other  second¬ 
ary  addressing  with  the  HP- 8 5  computer  requires  the  custom 
HP-IB  command  "SEND"  which  is  described  in  detail  in  the 
HP-IB  I/O  programming  manual  [Ref.  3] . 

The  12  bit  data  word  (voltage  reading)  is  routed  to 
the  A/D's  edge  connector  so  that  the  A/D  can  be  used  in  con¬ 
junction  with  the  69790B  memory  card  to  make  measurements  at 
rates  that  exceed  the  throughput  of  the  controller.  Up  to 
33000  readings  per  second  can  be  transferred  from  the  A/D 
converter  to  the  memory  card  via  the  edge  connectors  on  each 
card.  Each  reading  can  be  initiated  by  program  commands  or 
by  an  External  Trigger  (EXT)  signal. 

Either  the  69736A  Timer  Pacer  or  the  69735A  Pulse 
Train  Output  card  can  be  used  to  externally  trigger  the  A/D 
converter.  The  69736A  Timer/Pacer  is  programmed  in  a  con¬ 
venient  seconds-milliseconds-microseconds  format  for  pro¬ 
grammable  delay  or  continuous  triggering. 

Figure  3  shows  the  functional  block  diagram  of  the 
69751A  Analog  to  Digital  converter  card.  Various  subad¬ 
dresses  can  also  be  written  using  the  software  described  in 
section  5  of  the  User's  Guide  [Ref.  2]. 

3 .  Memory  Card  I69790B 

The  69790B  Memory  Card  is  a  FIFO  (first  in  first  out) 
buffer  that  can  handle  a  continuous  stream  of  data  that  is 
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being  transferred  from  the  outside  world  to  the  controller 
or  from  the  controller  to  the  outside  world.  In  another 
operating  mode,  the  recirculating  input  mode,  it  can  act  as 
a  circular  buffer  that  continuously  stores  data  coming  in 
from  the  outside  world  until  it  is  programmed  to  stop  ac¬ 
cepting  new  data.  At  this  time  it  contains  the  last  'n' 
words  of  data  received,  'n'  being  the  capacity  of  the  card. 

One  more  operating  mode,  the  recirculating  output  mode,  al¬ 
lows  all  or  part  of  the  memory  contents  previously  acquired 
to  be  read  repeatedly  by  an  external  device. 

Each  69790B  memory  card  consists  of  two  plug-in 
cards  that  occupy  two  adjacent  slots  in  the  mainframe.  The 
card  designated  card  #1  occupies  the  left-hand  slot  and  com¬ 
municates  with  card  #2  through  a  cable. 

A  memory  card  assembly  can  communicate  bidirection¬ 
ally  with  the  controller  and  with  an  external  device  con¬ 
nected  to  the  edge  connector  on  card  #1  of  the  pair.  External 
data  transfers  through  16  input  or  16  ouizput  lines  and  are 
controlled  by  means  of  three  input  or  three  output  handshake 
lines.  Special  purpose  Memory  Input  (MI)  and  Memory  Output 
(MO)  instructions  allow  data  transfers  to  be  made  between 
the  controller  and  a  memory  card  assembly  without  storing 
the  data  in  mainframe  memory. 

In  a  system  that  contains  more  than  one  memory  card 
assembly,  each  assembly  can  handle  data  transfers  to  or  from 
the  outside  world  simultaneously  with  external  data  transfers 
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being  made  by  the  other  assemblies.  Absolutely  no  inter¬ 
actions  occur  between  their  external  interfaces.  An  external 
interface  is  independent  not  only  of  the  other  external  in¬ 
terfaces,  but  also  of  any  other  card  or  mainframe  activity. 
Thus,  within  a  single  assembly,  an  external  device  can  write 
to  a  memory  card  while  the  controller  is  reading  from  another 
part  of  the  same  card,  or  an  external  device  can  read  from  a 
memory  card  while  the  controller  is  writing  to  another  part 
of  the  same  card.  Through  the  use  of  controller  interrupts 
or  another  means  of  program  control,  data  can  be  loaded  into 
memory  and  read  from  memory  on  a  continuous  basis  in  either 
an  input  or  an  output  operation. 

Figure  4  is  a  functional  block  diagram  of  the  69790B 


memory  card. 


Figure  4.  Mu ltiprogr amme r — Memory  Board  [Ref.  2] 


III.  PARTICLE  SIZE  DETERMINATION  BY  LIGHT  SCATTERING 
OF  A  HELIUM  NEON  LASER  BEAM 

The  purpose  of  this  experiment  is  to  determine  the  volume 
to  surface  mean  diameters  (D^)  particulate  matter  in  the 
exhaust  plume  and  combustion  chamber  of  a  solid  propellant 
rocket  motor. 

The  equipment  test  set  up  is  shown  in  Figure  5.  Refer¬ 
ence  4  provides  extensive  information  on  the  physical  princi¬ 
ples  involved  and  a  description  of  the  optics  system. 

The  helium  neon  laser  is  remotely  controlled  under  soft¬ 
ware  control  and  radiates  a  beam  at  a  wavelength  of  0.6328 
micrometers.  The  beam  is  then  expanded  and  columated.  It 
is  then  passed  through  a  limiting  orifice  to  reduce  the  beam 
diameter  to  2.2  millimeters. 

The  beam  is  subsequently  passed  into  a  beam  splitter 
where  one  beam  is  turned  through  90°  and  routed  through  a 
window  in  the  rocket  motor  casing.  The  other  beam  is  aligned 
directly  in  the  exhaust  stream  of  the  rocket  motor  approxi¬ 
mately  two  inches  behind  the  nozzle  throat. 

Each  beam  is  passed  through  a  focusing  lens  and  directed 
through  a  narrow  pass  filter  to  the  edge  of  a  linear  array 
of  photodiodes. 

Two  EG&G  RETICON  G  Series  Solid  State  Line  Scanners  are 
used,  each  containing  1024  photodiode  elements.  The  Reticon 


RL1024G  series  photodiode  array  is  optimized  for  solid  state 
image  sensor  applications.  It  contains  1024  diodes  arranged 
in  a  linear  array  on  25  micron  centers.  Diode  #1  is  normally 
set  1.54  millimeters  off  of  the  laser  beam  centerline  to  re¬ 
duce  the  amount  of  transmitted  light  present,  and  to  allow 
for  increased  sensitivity  to  the  diffracted  light  pattern. 

The  G  series  array  is  connected  through  an  RC-100/106 
sample  and  hold  "Boxcar"  circuit  which  provides  its  own 
clock  oscillator  for  self  scanning,  although  provisions  for 
externally  pacing  the  array  with  an  external  clock  signal  are 
installed. 

Photodiode  array  #1  (Fig.  5)  is  operated  in  the  internal 
clock/ internal  start  mode  and  is  completely  independent.  The 
internal  clock  oscillator  is  set  to  approximately  30  KHz.  To 
prevent  integrated  "dark  current"  from  making  a  significant 
contribution  to  the  output  charge  the  total  time  between 
start  pulses  should  be  kept  less  than  about  40  milliseconds. 

The  countdown  circuit  on  the  array  board  counts  ’n' 
clock  pulses  to  generate  a  start  pulse  in  the  internal  start 
mode.  The  countdown  switches  are  initially  set  according  to 
the  formula : 

n  ■  H  •  £c 

where  t^  *  time  between  start  pulses  in  seconds 
fc  »  clock  frequency  in  Hz 
n  *  switch  setting 
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e.g.  n  -  36(10"3)  •  30(103) 
n  -  1080 

A  countdown  switch  setting  of  1080  to  1  will  provide  a  time 
between  start  pulses  of  36  milliseconds  at  a  clock  frequency 
of  30  KHz,  and  excessive  "dark  current"  error  build-up  is 
avoided. 

Photodiode  array  #2  monitors  the  second  laser  beam  which 
emerges  from  a  window  in  the  rocket  motor  body.  Array  #2 
operates  in  the  external  clock/external  start  mode  and  is  a 
slave  to  photodiode  array  #1.  Therefore,  the  countdown  cir¬ 
cuits  in  array  #2  are  not  used.  Array  #2  has  only  one  out¬ 
put,  the  video  "boxcar”  signal  representing  the  photodiode 
voltage  output  which  is  routed  directly  to  A/D  converter  #2 
in  the  6942A  multiprogrammer. 

The  boxcar  voltage  outputs  from  the  photodiode  arrays 
are  routed  through  analog  low  pass  filters  which  pass  all 
signal  components  less  than  3000  Hz.  20  db  gain  is  attained 
in  the  low  pass  filter  amplifier.  The  filtered  outputs  are 
routed  to  the  voltage  input  terminals  (+W,  -Y)  of  the  high 
speed  analog  to  digital  converters  located  in  the  6942A 
multiprogramraer .  When  a  series  of  readings  are  desired  the 
HP-85  programs  the  69736A  timer  to  output  a  square  wave  output 
at  approximately  30  KHz  when  externally  triggered  by  the 
blanking  pulse  output  from  photodiode  array  #1. 

The  square  wave  output  from  the  timer  pacer  card  is  di¬ 
rected  to  the  external  trigger  inputs  of  both  69751A  high 
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speed  (A/D)  cards  which  begin  making  1024  conversions  of  the 
analog  voltage  present  at  terminals  (+W,  and  -Y) .  The  12  bit 
digital  outputs  from  terminals  A  through  N  are  hard  wired  to 
the  corresponding  edge  connector  pins  of  the  69790B,  4K  memo¬ 
ry  cards  associated  with  each  (A/D)  card.  The  memory  is  pro¬ 
grammed  to  accept  a  12  bit  format  and  to  store  1024  readings 
in  a  first  in  first  out  (FIFO)  input  mode.  The  memory  card 
is  then  armed  to  interrupt  the  computer  after  1024  readings 
have  been  stored  on  the  card.  The  (A/D)  card  performs  a  con¬ 
version  of  the  analog  voltage  at  its  input  terminals  (W,  Y) 
every  time  it  observes  a  trigger  from  the  timer  pacer  card  at 
its  external  trigger  terminals  (m)  at  a  rate  of  30  KHz.  The 
(A/D)  card  signals  the  completion  of  each  reading  by  output¬ 
ting  an  end  of  conversion  (EOC)  and  its  complement  (EOC) . 

The  E6C  is  connected  to  the  data  available  (DAV)  terminal  of 
each  69790B  memory  board. 

The  procedure  employed  for  an  experiment  is  outlined  be¬ 
low.  Figure  6  presents  a  flow  diagram  of  the  system.  After 
applying  power,  the  HP-85  is  loaded  with  a  program  disc  in 
drive  #0  and  a  data  disc  in  drive  #1.  The  program  named 
"HELP"  is  loaded  and  run  to  initialize  variables  and  provide 
documentation.  At  the  completion  of  "HELP"  the  next  program 
is  automatically  loaded  and  executed  by  a  process  called 
chaining,  which  allows  a  program  of  virtually  unlimited  length 
to  be  run  in  the  limited  memory  space  of  32K  of  the  HP-85 
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Figure  6 
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computer.  "ACQDTA”  now  performs  the  process  control,  system 
timing,  and  data  acquisition  tasks. 

First,  the  laser  is  energized,  data  files  are  created 
with  user  specified  file  names,  and  initial  readings  are 
taken  from  both  photodiode  arrays  by  software  command.  Each 
array  simultaneously  sends  its  output  to  its  69751A  (A/D) 
converter  in  the  6942A  multiprogrammer.  Each  (A/D)  card  is 
externally  triggered  every  33  microseconds  by  the  69736A 
Timer/Pacer  card  and  converts  the  analog  voltage  present 
at  terminals  (+W,  -Y) .  The  resulting  12  bit  digital  out¬ 
put  is  sent  to  a  69790B  (4K)  memory  card.  Each  memory  card 
is  programmed  to  store  sequentially  1024  readings  and  to 
interrupt  the  computer  after  1024  readings  have  been  stored. 
At  this  time  the  numerical  data  are  transferred  via  the 
HP- IB  to  the  controller,  first  array  #1  then  array  #2.  This 
raw  data  represents  the  "no  particle"  or  "zero  scattering" 
condition. 

The  "no  particle"  numerical  arrays  are  then  stored  as 
a  data  file  on  the  disc  to  preclude  loss  of  the  data. 

The  computer  then  performs  a  continuity  check  of  the 
firing  circuit  using  the  3456A  in  a  resistance  measurement 
mode.  After  a  successful  continuity  check,  the  3456A  is  pro¬ 
grammed  to  record  DC  voltage.  It  is  connected  to  the  firing 
switch  while  displaying  (on  the  HP-85) ,  "STANDING  BY  FOR 
IGNITION",  and  energizes  the  nitrogen  purge  for  the  motor 
windows.  When  the  operator  closes  the  ignition  switch  a 


voltage  greater  than  10  volts  is  sensed  by  the  3456A  and 
causes  the  computer  to  energize  the  drive  motor  of  the  visi¬ 
corder,  which  records  the  chamber  pressure-time  trace. 

Once  the  computer  senses  the  ignition  sequence  it  con¬ 
nects  the  3456A  to  the  pressure  transducer  in  parallel  with 
the  visicorder  trace. 

When  the  combustion  pressure  rises  to  a  user  specified 
threshold  pressure  in  (PSI)  the  computer  delays  or  waits  a 
user  specified  number  of  seconds,  and  triggers  the  multipro¬ 
grammer  to  acquire  a  new  data  set  from  the  photodiode  arrays. 
When  the  data  are  taken  a  pulse  output  voltage  (TTL)  is  sent 
to  another  visicorder  channel  as  an  event  marker  signal  to 
indicate  data  acquisition  time  on  the  pressure-time  trace. 
Readings  taken  during  the  rocket  motor  burn  time  will  be 
different  than  pretest  readings  since  particulate  matter  is 
present  in  the  gas  stream.  This  data  is  transferred  to  the 
computer  and  then  to  a  data  file  on  the  data  disc  as  before. 

After  5  seconds  the  visicorder  is  disengaged,  the  window 
nitrogen  purge  turned  off,  and  the  laser  is  deenergized.  The 
next  program  is  chained. 

"CMBINE"  is  automatically  loaded  and  run.  This  program 
executes  disc  operations  which  read  the  "no  particle"  data 
and  the  "particle"  data  into  the  computer.  It  then  separates 
this  data  into  photodiode  array  #1  data  and  photodiode  array 
#2  data  with  the  same  file  names  as  previously  specified  by 
the  user.  Each  file  contains  1024  records,  each  containing 
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two  numbers,  a  "no  particle"  and  a  "particle"  voltage.  The 
computer  prints  a  reminder  of  which  data  is  in  which  data 
file  and  then  chains  the  next  program  segment. 

"PLTDTA"  plots  the  data  from  a  data  file  specified  by 
the  user,  plotting  diode  voltage  output  versus  diode  number 
(Fig.  7) .  The  operator  cam  elect  to  plot  the  other  data 
file  or  load  and  run  the  data  reduction  prograun. 

The  operator  is  then  prompted  to  load  and  run  "REDUCE"  to 
begin  the  reduction  of  the  raw  data.  Figure  8  shows  a  flow 
chart  for  the  data  reduction  process.  "Reduce"  functions  to 
compute  a  normalized  intensity  vector  and  an  angle  theta 
vector.  The  theta  vector  is  computed  based  on  diode  number 
(1  to  1024)  ,  lens  focal  length,  and  the  index  of  refraction 
of  the  media  containing  the  particles  by  the  formula: 

K1  +  n ( .025) 

9  =  - - - -  radians 

f  m 

where  K.  *  the  distance  from  the  centerline  of  the  laser  beam 
to  the  centerline  of  diode  #1  (mm) 

n  =  diode  number  (1  to  1024) 

F  *  lens  focal  length  (mm) 

m  *  index  of  refraction 
air  *  1.0 
water  =*  1.35 

.025  *  distance  between  diode  cent...  s  (mm) 

Upon  examination  of  the  "PLTDTA"  results,  the  location 
of  the  peak  intensities  are  printed.  The  user  is  then  prompt¬ 
ed  for  the  diode  locations  where  the  slope  is  to  be  determined 
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Figure  8.  Light  Scattering  Program  Flow  Chart — Part  Two 


and  then  projected  back  to  obtain  a  centerline  value  for  the 
scattered  laser  power  using  a  linear  regression  technique. 

First,  the  "zero  scattering"  data  is  subtracted  from  the 
"with  particle"  data.  The  resulting  1024  element  vector  is 
then  regressed  against  diode  number  to  obtain  the  intensity- 
axis  intercept  representing  the  centerline  value  by  the  fol¬ 
lowing  formula: 

nA  +  BX1  =  Y1 
xxA  +  BX2  =  Y2 


where  i  *  1  to  1024 

IqNOrm  “  normalized  intensities 
VQp(i)  ■  voltage  with  particles  (vector) 
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VeNp(i)  *  voltage  without  particles  (vector) 

(vCLp  “  vCLfjp)  *  A  (previous  calculation) 

The  normalized  intensity  is  then  plotted  on  the  7225B 
plotter  versus  the  theta  vector  in  radians  (Fig.  9) .  This 
data  is  then  stored  on  disc  for  future  use  in  a  user  speci¬ 
fied  file.  The  operator  is  then  asked  if  this  run  is  to  be 
a  calibration  curve;  if  yes,  then  a  curve  fitting  program 
"CALIB"  is  automatically  loaded  and  executed. 

"CALIB"  performs  a  tenth  order  polynomial  curve  fit  to 
the  data  set  specified  by  the  user  in  order  to  generate  a 
smooth  curve  for  future  comparison  purposes  with  both  the 
universal  (i.e.  theoretical  [Ref. 4])  curve  ("UNICUR")  and 
other  raw  data  runs.  The  user  first  enters  the  "D  bar" 
value,  i.e.  the  average  diameter  of  the  particles  being  used 
for  the  calibration. 

A  theta  bar  vector  is  then  computed  from  the  theta  vector 
as  follows: 


where  9^  =  theta  vector  (1-1024)  (radians) 

D32  *  D  bar  (microns)  mean  volume  to  surface  diameter 
X  =  .6328  (microns)  wavelength 
Normalized  intensity  is  then  plotted  versus  the  non- 
dimensional  theta  bar  (9)  as  in  figure  10.  The  user  is  then 
asked  to  enter  the  minimum  and  the  maximum  theta  bar  of 
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THETA  (RADIANS) 

REDUCE  Output  Normalized  Intensity  vs  Theta 


interest.  The  data  is  regressed  to  form  a  smooth  curve  fit 
with  a  polynomial  of  up  to  tenth  order  as  specified  by  the 
user  and  taking  the  form: 

IN  -  b„  .  bl6  -  V*  ♦  b383  ♦  .  .  .  +  bme» 

where  the  coefficients  are  obtained  by  solving  the  system  of 
m+1  equations 


V9 i”  +  bil9"*1  +  •  •  .  +  bje2”  *  Ili9m 

Once  the  polynomial  coefficients  have  been  determined,  a 
smooth  100  element  curve  is  computed,  plotted  and  stored  on 
the  disc  in  a  user  specified  file. 

The  smooth  curve  is  plotted  on  top  of  the  intensity  vs 
theta  bar  plot  (Fig.  10)  to  show  the  correlation,  which  is 
routinely  excellent.  If  desired,  this  data  file  containing 
smooth  curve  data  can  be  saved  permanently  on  the  floppy  disc. 

"REDUC2"  is  automatically  loaded  and  executed  upon  com¬ 
pletion  of  "CALIB".  The  function  of  "REDUC2"  is  to  plot  the 
universal  theoretical  curve  (datafile  "UNICUR" )  on  a  three 
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cycle  semi- logarithmic  scale  (Fig.  11) .  After  the  theoretical 
curve  is  plotted  the  user  may  elect  to  plot  any  other  calibra¬ 
tion  curve  for  comparison  purposes  (Fig.  11) . 

At  the  conclusion  of  NREDUCE",  if  the  operator  had  not 
designated  a  run  as  a  calibration  run,  the  next  data  reduc¬ 
tion  program  is  chained. 

"REDUC3"  functions  to  compute  from  the  normalized 
intensity  and  theta  data  by  comparing  the  data  to  specified 
reference  curves  (UNICUR,  etc.)  at  15  preset  intensity  levels. 
An  average  is  computed  across  this  range  of  intensity 
levels.  REDUC3  has  tabular  output  only. 

REDUC4  can  be  run,  if  desired,  to  plot  a  data  run  on  a 
semi- log  scale  with  a  user  defined  D32  to  compare  with  the 
Universal  curve. 
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VITIATED  AIR  HEATER 


A  schematic  of  the  vitiated  air  heater  system  at  the  com¬ 
bustion  laboratory  is  shown  in  figure  #12.  It  was  designed 
to  supply  air  at  any  desired  mass  flow  rate  within  the  system 
limitations  (approximately  4  lbm/sec)  and  at  any  temperature 
up  to  approximately  1500°R.  The  heater  output  air  is  routed 
to  the  input  sections  of  the  Solid  Fuel  Ramjet  or  the  Turbo¬ 
jet  combustor. 

The  vitiated  air  heater  is  fueled  by  gaseous  ethylene  and 
bottled  oxygen  replaces  that  consumed  in  the  combustion  pro¬ 
cess.  Large  high  pressure  air  tanks  act  as  the  primary  air 
supply  and  the  mass  flow  rate  is  determined  using  a  sonically 
choked  converging  nozzle. 

The  data  acquisition  system  was  designed  to  monitor  and 
control  the  operation  of  the  air  heater  in  order  to  provide 
the  desired  air  flow  rate  and  temperature  and  to  provide 
built-in  system  safety. 

Two  Hastings  precision  flowmeters  and  two  Hastings  flow 
controllers  were  installed  to  accurately  meter  and  control 
the  oxygen  and  ethylene  flow  rates  to  achieve  the  desired 
fuel-air  ratio  and  oxygen  make-up  flow  rate.  Figure  #13 
shows  the  interconnection  of  the  flowmeters  and  controllers. 

A  description  of  the  system  operation  is  presented  below. 


The  'HEATER'  program  disc  is  loaded  into  the  HP-85  data 


acquisition  system.  The  "SET  UP”  option  is  selected  to  ini¬ 
tialize  variables  and  perform  preliminary  calculations  to 
determine  initial  manual  control  settings.  Figure  14  depicts 
the  flow  chart  for  the  preliminary  portion  of  the  ’HEATER' 
program. 

The  operator  inputs  the  desired  mass  flow  rate  of  air  in 
lbs  mass/ sec.  Other  inputs  are  high  pressure  air  tank  pres¬ 
sure  (0-3500  PS1G) ,  the  estimated  cold  air  temperature  (°F) , 
desired  hot  air  temperature  (°R) ,  and  ambient  pressure  (inches 


The  required  fuel/air  ratio  is  then  computed  for  ethylene 
and  air  from  a  regression  of  data  obtained  from  the  "PEPCODE" 
chemical  equilibrium/adiabatic  combustion  program  [Ref.  5] . 

„  _  TH0T  “  550  *FUEL 
78000  " 


Once  F  is  obtained,  the  ethylene  and  oxygen  flow  rates 


are  given  by: 


m„  „  =  F  •  m_  (#  mass/sec) 

c2h4  °2 


For  the  complete  combustion  of  ethylene  with  oxygen  in  air: 

C2H4  +  302  -  2C02  +  2H20 

Thus,  the  desired  mass  flow  ratio  of  oxygen  to  ethylene  is 

”o2  "  3-429  •  "c2h4 


'  i 

•  * 


Figure  14.  HEATER  Flow  Chart — Part  One 


The  volume  flow  rates  of  oxygen  and  ethylene  in  SCFM  are: 


Q.  =  734.89  •  mn 
°2  u2 


V4  -  339'87  •  ”c2H4 


thus 


°°2  '  3V< 


where  Q  is  the  volume  flow  rate.  The  oxygen  and  ethylene 
flowmeters  produce  0-5  Volt  DC  signals  linear  with  volume 
flow  rate  (0-20  SCFM  for  oxygen  and  0-6.5  SCFM  for  ethylene). 
The  required  oxygen  and  ethylene  set  voltages  are  computed 
from 

V_  =  183723  •  mA  (lbm/sec)  (millivolts) 
u2  u2 

VA  „  »  646054  •  iru  u  (lbm/sec)  (millivolts) 

VQ  is  generated  by  the  dual  channel  digital  to  analog 
(D/A)  converter  card  in  the  3497A  and  sent  to  the  REMOTE 
terminals  of  the  oxygen  flow  controller  (Fig.  13) .  This  flow 
controller  is  set  to  the  'BIAS'  mode  of  operation  with  the 
Local  Set  Point  (LSP)  control  set  to  zero. 

In  this  position  (Fig.  15a)  the  input  stages  of  the  flow 
controller  sum  the  three  voltages  and  generate  an  error 
signal : 


ERROR  =  INPUT  -  LSP-  REMOTE 


This  signal  drives  the  flow  control  until  the  error  signal  is 
equal  to  zero,  which  occurs  when  the  desired  oxygen  flow  rate 


is  achieved.  With  the  LSP  set  to  000,  the  oxygen  set  voltage 
from  the  D/A  converter  acts  as  a  remote  set  point  to  control 
the  desired  flow  rate.  The  upper  section  of  figure  13  shows 
the  interconnections  of  the  second  flow  controller,  the  ethyl¬ 
ene  flowmeter,  and  the  oxygen  flowmeter  in  the  RATIO  mode  of 
operation.  In  this  mode,  the  flow  controller  acts  to  provide 
the  proper  ratio  of  ethylene  to  oxygen  under  all  conditions 
by  controlling  the  ethylene  volume  flow  rate  to  maintain  a 
1:3  ratio.  Thus,  the  computer  controls  the  oxygen  flow  rate 
which  in  turn  controls  the  ethylene  flow  rate. 

Figure  15b  shows  the  flow  controller  input  circuit  in  the 
1-5  RATIO  mode.  In  this  case,  the  REMOTE  (OXYGEN)  signal  is 
equal  to  the  product  of  the  Local  Set  Point  (LSP)  and  the 
INPUT  (ETHYLENE)  signal. 

LSP  •  INPUT  *  REMOTE 

However,  the  input  (ethylene)  and  remote  (oxygen)  flow 
capacities  are  different.  Dividing  by  the  full  scale  flow 
capacities  in  SCFM  normalizes  the  equation. 

LSP  •  1-575“/  =  \  15  / 

REMOTE  =  Oxygen  volume  flow  rate 
INPUT  =  Ethylene  volume  flow  rate 

Then, 

INPUT  _  l  6 . 5 \/  1  \  _  1 
REMOTE  \  20  /\LSP/  "  T 


Thus,  LSP  =  0.975. 


For  example,  if  the  desired  oxygen  flow  rate  is  3.5  SCFM 
or  0.00476  lbm/sec,  then  the  OXYGEN  set  voltage  is 

»0j  -  3.5  SCFM  foV|^  -  0.875  Volts 

The  computer  sends  0.875  volts  to  the  remote  terminals  of 
the  oxygen  flow  controller.  The  RATIO  controller  responds  by 
driving  the  flow  control  valve  until  the  ethylene  flowmeter 
reads: 

INPUT  ♦  LSP  =  REMOTE 
INPUT  =  =  0.8974  VOLTS 

This  corresponds  to  an  ethylene  volume  flow  rate  of 

Qc2h4  =  °‘8974  Volts 

Q_  „  =1.166  SCFM 

C2H4 

which  is  1/3  of  the  oxygen  volume  flow  rate. 

This  completes  the  preliminary  calculations  and  the  sys¬ 
tem  is  now  turned  on  to  provide  air  flow  without  any  heat 
addition.  Thus,  the  operator  selects  the  'FLOW'  option  (Fig. 
14)  and  the  computer  measures  ambient  pressure  then  opens  the 
air  valve,  allowing  air  to  flow  to  the  heater  while  other 
components  are  maintained  in  the  off  position.  High  pressure 
air  passes  through  a  choked  flow  nozzle,  where  PTl  and  T^ 
are  measured.  Based  on  these  measurements  the  actual  mass 
flow  rate  is  computed  from: 
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(lbm/sec) 


CD  PTl  A*  F1 


where 


F 


1 


Y 


A* 

d* 

CD 

R 


.235  (in)  typical 
0.97,  discharge  coef. 


ft-lbF 

lbm4R 


T<j<1  *  measured  temperature,  °R 


The  actual  mass  flow  rate  of  air  is  displayed  and  an  up¬ 
dated  oxygen  set  voltage  is  computed  and  set.  The  required 
pressure  correction  is  computed  to  meet  the  desired 

mass  flow  rate. 


The  operator  may  then  accept  the  flow  conditions  and 
proceed  to  another  program  segment.  If  he  accepts  these 
flow  conditions,  the  actual  mass  flow  rate  of  air  is  set 


equal  to  the  desired  mass  flow  rate.  If  he  is  not  satisfied 
with  existing  flow  conditions,  is  adjusted  manually  by 


56 


applying  the  recommended  AP^ •  The  'FLOW'  option  is  auto¬ 
matically  executed  again  until  the  desired  accuracy  is 
attained. 

Once  the  desired  mass  flow  rate  of  air  is  obtained,  the 
'START'  option  is  executed.  Figure  16  shows  a  flow  chart  of 
this  program  segment.  This  routine  first  closes  all  valves 
and  measures  PTl  which  represents  ambient  pressure  under  no 
flow  conditions.  The  proper  oxygen  set  voltage  is  set  ac¬ 
cording  to  the  latest  flow  calculations.  The  computer  then 
turns  on  the  high  pressure  air  flow,  energizes  the  igniter 
system,  and  turns  on  the  ethylene  and  oxygen  flow  to  the 
heater  section.  The  output  hot  air  temperature  is  monitored 
by  a  thermocouple  and  measured  by  the  computer  at  0.5  second 
intervals  during  the  start  up  sequence,  beginning  3  seconds 
after  ignition.  The  igniter  is  turned  off  after  5  seconds  of 
operation  unless  secured  earlier  by  manual  control. 

During  heater  start  up  a  typical  temperature-time  trace 
would  be  as  shown  in  figure  17. 

Temperature  is  closely  monitored  during  the  start  se¬ 
quence  and  igniter  shutdown  is  executed: 

1)  when  THqT  air  reaches  700°R  (if  not  previously  secured 
by  the  5  second  timer) 

2)  if  TgoT  AIR  falls  by  50 °R  in  any  5  second  interval  then 
HEATER  shutdown  is  called  to  close  the  OXYGEN  and  ETHYL¬ 
ENE  valves,  close  igniter  valves,  set  O2  set  voltage  *  0 
and  close  the  high  pressure  air  valve. 

When  Thqt  air  reaches  TDESIRED  -  100°  the  Heater  Mainten¬ 
ance  routine  is  called,  terminating  the  start  sequence.  The 
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Figtire  16.  HEATER  Flow  Chart — Part  Two 


Heater  Maintenance  Program  is  assigned  to  timer  #3  and  set  to 
automatically  execute  every  20  seconds,  while  another  project 
program  segment  is  run. 

Figure  18  shows  the  flow  chart  for  the  Heater  Maintenance 
Program  segment  (HTRMNT) ,  which  is  assigned  to  Timer  #3  to 
interrupt  the  computer  at  20  second  intervals  to  control  the 
heater  output  temperature.  On  interruption,  the  computer 
measures  the  output  air  temperature  and  compares  it  to  the 
desired  air  temperature.  If  the  difference  is  greater  than 
200°R  the  heater  shutdown  program  is  called  and  secures  the 
heater  entirely.  If  the  difference  is  less  than  100°R  the 
computer  records  the  air  temperature  and  returns  to  whatever 
routine  was  operating  prior  to  the  interrupt.  However,  if 
the  difference  between  actual  temperature  and  desired  tempera¬ 
ture  is  over  100°  and  less  than  200 °R  the  sign  of  the  error 
is  determined  and  the  fuel  air  ratio  is  adjusted  to  compen¬ 
sate  for  the  temperature  deviation.  Oxygen  and  ethylene  flow 
rates  are  updated  and  a  new  oxygen  set  voltage  is  sent  to  the 
oxygen  flow  controller  to  effect  a  change  in  the  output  tem¬ 
perature.  The  output  temperature  is  then  recorded  and  a  re¬ 
turn  from  interrupt  is  executed. 
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V.  SOLID  FUEL  RAMJET 


The  Solid  Fuel  Ramjet  (SFRJ)  test  facility  at  the  combus¬ 
tion  laboratory  has  been  constructed  to  facilitate  the  mea¬ 
surement  of  performance  as  affected  by  design  and  operating 
variables.  Figure  19  presents  a  schematic  of  the  SFRJ  motor 
which  is  mounted  on  a  thrust  stand.  Figure  20  presents  a 
schematic  of  the  air  supply  system.  Air  flow  is  directed  to 
the  head-end  and  sometimes  also  to  the  aft  mixing  chamber 
(bypass-air)  of  the  combustor,  in  general,  a  test  sequence 
is  as  follows:  air  flow  rate(s)  are  turned  on,  the  igniter 
is  actuated,  a  small  amount  of  ethylene  is  mixed  with  the 
head-end  air  flow,  the  motor  ignites,  the  igniter  and  ethyl¬ 
ene  flows  are  terminated,  the  motor  runs  for  10-45  seconds 
duration,  air  flow  is  terminated  and  nitrogen  gas  is  flowed 
through  the  fuel  grain  for  five  seconc-  to  ensure  extinguish¬ 
ment  of  the  flame.  During  the  test  it  is  important  to  mea¬ 
sure  air  flow  rates,  inlet  air  temperature,  head-end,  and 
mixing  chamber  pressures.  Postfire  data  collected  include 
fuel  grain  weight  loss  which  is  needed  to  calculate  the 
average  fuel  flow  rate  and  average  fuel  regression  rate. 

Temperature  rise  combustion  efficiency  for  the  test  is 
determined  from 


AT  Tv*.  -  Tt  . 

^th  '■air 


j -i  *mrn  ■' 


Figure  20.  Solid  Fuel  Ramjet  Functional  Block  Diagram 


where  Tta^r  is  the  measured  inlet  air  stagnation  temperature, 
Ttt^  is  the  theoretical  combustion  temperature  determined 
from  adiabatic,  equilibrium  combustion  calculations  using 
the  measured  fuel-air  ratio,  inlet  air  temperature  and  aver¬ 
age  chamber  pressure.  TteXp  is  the  experimentally  determined 
combustion  temperature.  The  latter  is  calculated  from  con¬ 
servation  of  mass,  using  gas  properties  determined  from  the 
equilibrium  calculations.  Thus, 


Pc  =  average  chamber  stagnation  pressure 
CD  =  nozzle  discharge  coefficient 


Average  fuel  regression  rate  is  calculated  based  upon 
fuel  weight  loss  during  the  burn.  Thus, 


where  AW  *  fuel  weight  loss 
L  =  fuel  grain  length 


- - — - -  11 1  '1 

p  =  fuel  density 

=  final  average  fuel  internal  diameter 

t.  =  burn  time 
b 

J  =  average  fuel  regression  rate 

The  average  chamber  pressure  (Pc)  used  to  determine  Ttexp 
is  obtained  from  the  pressure  time  trace  recorded  during  the 
test.  A  typical  trace  is  shown  in  figure  21. 

The  test  control/data  acquisition/data  reduction  program 
that  performs  the  above  experiment  is  described  in  this 
section. 

Figure  22  depicts  the  flow  chart  for  the  Solid  Fuel  Ram¬ 
jet  ' SFRJ '  option,  which  is  assigned  to  special  function  key 
#3  on  the  HP-85.  This  option  is  executed  after  having  loaded 
and  run  the  'HEATER'  program  (figure  14) ,  completed  the 
'SETUP'  option,  and  executed  the  'FLOW'  option  to  set  up  the 
proper  flow  rates  of  air,  oxygen,  and  ethylene.  If  the 
vitiated  air  heater  is  to  be  used,  the  'START'  option  (fig¬ 
ure  16)  is  executed  to  bring  ;he  heater  up  to  stable  operat¬ 
ing  temperature.  Once  the  air  temperature  and  total  air  flow 
rate  (ifi^)  are  within  desired  limits  of  the  required  values, 
the  heater  maintenance  'HTRMNT'  (figure  18)  is  automatically 
set  to  interrupt  the  computer  at  10  second  intervals  to 
check  the  operation  of  the  heater. 

Once  the  'SFRJ'  key  is  pressed,  the  operation  of  the 
Solid  Fuel  Ramjet  is  initiated.  Both  air  flow  control  valves 
are  initially  opened  to  approximately  the  50%  position. 
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TIME  (SECONDS) 

#21  Typical  Pressure  vs.  Tine  Trace. 


Figure  22.  Solid  Fuel  Ramjet  Flow  Chart 
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Figure  22  (Continued) 
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Initially,  heater  output  air  is  routed  via  primary  and  bypass 
a^r  lines  to  air  dump  valves,  which  are  held  in  the  dump  po¬ 
sition,  venting  heater  air  to  atmosphere.  The  SFRJ  igniter 
and  the  ethylene  valves  are  also  maintained  off  at  this  time. 

The  primary  valve  position  is  then  adjusted  by  the  com¬ 
puter  until  the  primary  air  flow  rate  (ih^)  is  within  5%  of 
the  desired  value.  The  computer  programs  the  6942A,  TIMER/ 
PACER,  A/D  converter,  and  memory  cards  to  read  the  SFRJ  motor 
chamber  pressure  (Pc)  transducer  output  voltage  continuously 
for  60  seconds  after  being  externally  triggered.  The  reading 
rate  is  fixed  at  approximately  16  readings  per  second. 

When  the  6942  is  programmed  the  SFRJ  manual  ignition 
switch  is  enabled  by  the  computer  which  then  monitors  this 
switch  until  it  detects  closure.  Once  the  switch  is  closed 
and  sensed,  the  Timer /Pacer  card  is  triggered  and  the  com¬ 
puter  timer  is  set  to  zero.  Heater  air  is  then  directed 
into  the  SFRJ  motor,  the  igniter  is  turned  on,  and  the  ^2^4 
is  turned  on. 

The  computer  waits  a  preset  time  in  seconds  (TIGN)  before 
shutting  off  the  igniter  circuit  and  the  ^2^4  valve.  During 
the  start  sequence  the  3456A  digital  voltmeter  is  set  to 
monitor  chamber  pressure  (P^,)  .  The  chamber  pressure  is  mea¬ 
sured  and  compared  to  the  desired  steady  state  value.  If 
this  pressure  exceeds  .75  of  the  desired  value  by  the  end  of 
ignition  the  computer  proceeds  to  a  measurement  sequence  of 
several  pressures,  temperatures,  and  flow  rates  until  the 
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desired  burn  time  in  seconds  expires.  However,  if  the  pres¬ 
sure  is  less  than  .75  of  the  desired  value  the  program  di¬ 
rects  the  user  back  to  attempt  another  ignition  sequence. 


The  automatic  measurement  sequence  proceeds  from  the 
point  where  the  pressure  exceeds  .75  Pc  aesire(j  as  (refer 
to  figure  21): 

a)  P-T]/  tTi  are  measured 


b) 


m_ 


air 


pTx  A  Cp  Fl 
Q 


where  CD  =  .97 

A  =  ird2/4 
d  =  .235 


Y  =  1 . 4 


c) 

d) 


e) 

f) 

g) 


R  TTi 


P-r2 ,  TT2  measured  at  primary  air  orifice 
PT2  A2  CD  F1 


"P 


where  Q 


tT2 


Vm  is  set  to  obtain  desired  in 

mp  p 

•  •  • 

ItlBYPASS  "  mA  "'PRIMARY 

^BYPASS  =  0  no  bypass  selected 


THOT  AIR  ”easured 


h)  Tgp  measured 


i) 

j) 


SFRJ  head  pressure  and  temperature  measured  PH»  TH 

mp  TH  +  mgp  Tbp 


k)  Values  of  time,  maj.r,  T,  and  Tjjot  are  stored  in  an 
array  for  later  use  in  reduction  of  data. 

The  above  measurements  proceed  for  a  total  of  100  read¬ 
ings  at  a  rate  of  approximately  1.5  per  second. 

Both  the  SFRJ  and  the  heater  are  shut  down  automatically 
when  timer  #1  expires  at  (t_)  seconds  after  successful 

Q 

ignition. 

Once  the  measurement  sequence  completes  and  the  shutdown 
sequence  is  completed  the  computer  then  waits  until  the  multi¬ 
programmer  memory  board  signals  completion  of  its  storage 
operation. 

During  the  operation  of  the  SFRJ  option,  special  function 
key  #4  is  assigned  to  the  SFRJ  shutdown  routine.  This  pro¬ 
gram  is  automatically  executed  t^  seconds  after  ignition, 
however,  it  can  be  invoked  at  any  time  by  the  operator  and 
when  called  performs  the  following  steps: 

a)  Air  valve  to  dump  position 

b)  Nitrogen  purge  on  for  5  seconds 

c)  Heater  shutdown  sequence. 

After  all  timed  events  have  completed  the  computer  auto¬ 


matically  loads  and  executes  'SFRJ2',  which  is  the  data 


reduction  program  for  the  Solid  Fuel  Ramjet  experiment.  Fig¬ 
ure  23  shows  the  flowchart  of  the  data  reduction  program. 

A  large  quantity  of  measured  data  is  automatically  passed 
from  the  Heater  program  to  the  reduction  program.  The  com¬ 
puter  then  prompts  the  user  for  other  required  input  data, 
such  as: 

a)  fuel  type  HTPB/PMM 

b)  final  grain  weight  (grams) 

c)  initial  grain  weight  (grams) 

d)  initial  grain  internal  diameter  (cm) 

e)  final  aft-end  internal  diameter  (cm) 

f)  grain  length  (cm) 

2 

g)  area  of  the  exhaust  nozzle  throat  (cm  ) 

With  this  data  entered  and  stored,  the  computer  then  reads 
in  the  chamber  pressure  (Pc)  array,  which  is  a  1000  element 
vector  of  transducer  voltages.  These  voltages  are  subsequent¬ 
ly  converted  to  pressure  (PSIA)  by  the  formula 

PC(I)  -  K:(V(I)  -  V0 )  +  B9 

where  *  79.434 

V ( I )  *  transducer  voltage 

V0  *  initial  transducer  voltage  at  zero  psig 
B9  ■  atmospheric  pressure 

The  computer  then  locates  the  maximum  pressure  reading 
and  its  location  in  the  array.  The  computer  then  searches 
both  forward  and  backward  in  time  from  the  location  of  the 


INPUT  DATA 

1)  Enter  Init.  Grain  Wti 

2)  Enter  Final  Grain  WtJ 

3)  Enter  Grain  Langtb 

4)  Enter  Grain  Diaaeter 

5)  Enter  Final  Aft  Diask] 

6)  Enter  Fuel  Type - 


1)  Read  in  Pc(t) 

2)  Convert  P(t)  to  PSIG 


nr 

Find 

P-(max) 

2) 

Find 

Time  (Pc  max) 

3) 

Find 

Ts«Tbne(.  2*PC  max) 

4) 

Find 

Tp-TlnieC2^c  max) 

5) 

Tb  “ 

Tp  -  TS 

Int| 


GRATE 


J  Pc(t)dt 


P  m 

*C 


“  1  fF  . 

mA  ■  —  J  m.dt 
A  tb  tsJ  A 


AVE 


r?i 

i«r  +  di 


dp 

rave  i  ttLo 

7  _  ^Fave  “ 

ave  It7 

AWt  b 

°F  "TT 

f  -  mp/mA 

<t>  -  f/fST0ICH 


Figure  23.  Solid  Fuel  Ramjet  Reduction  Flow  Chart 
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PEPCODE  REGRESSIONS 
From  T"a.  and  Pq 
Calculate  Ttth,  Rc  and  Y 


Figure  23  (Continued) 


maximum  pressure  until  it  finds  the  time  where  pressure  has 
dropped  to  0.2  times  the  maximum  pressure  and  records  these 
locations  as  t  t  and  fcfiniah*  The  effective  final  burn 
time  is  tfiF  »  tp  -  tg.  The  computer  then  performs  a  numerical 
integration  of  the  pressure  time  trace  from  tg  to  tp  and 
divides  the  result  by  tBp  to  give  Pc,  the  average  pressure. 


Now  that  tg  and  tp  are  known,  the  computer  can  locate 
these  times  in  the  time  column  of  the  data  array  (100x3 
matrix)  and  note  the  row  numbers  of  I  and  I  .  The  average 
air  temperature  T_„_  and  flow  rate  rft  . _  are  then  computed 

ave  all 

between  tg  and  tp  by: 


*F 

E  *A 

S  .  = 

air  (IF-IS) 
JF 

E  T 

i=IS 

Tav  =  (IF-IS) 


Data  reduction  can  now  proceed  as  follows: 


AW. 


Initial  Wfc  -  Final  Wfc 


/ 4  AW.  ' 

/ - r  +  d." 

\  1 1  op  L  i 
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1 


where  pp  *  Fuel  propellant  density  (gms/cm3) 
L  *  length  (cm) 

*  initial  diameter  (cm) . 

The  average  burn  rate  is  then 


r 


ave 


The  equivalent  fuel  flow  rate 


“P 


Thus,  fuel/air  and  equivalence  ratios  are 


f 


=r 


5 

*A 


and 


y  -  j 

STOICH 

Then  with  T,  ,  and  the  type  of  fuel  burned,  the  theoretical 
values  of  Ttt^/  R^,  and  are  computed  from  a  polynomial 
regression  obtained  from  the  PEPCODE  [Ref.  5]  computer  pro¬ 
gram  at  the  Naval  Postgraduate  School  computer  center.  The 
effect  of  Pc  on  Ttth  is  negligible  for  nominal  variations 
(40-60psia) .  However,  Pc  effects  can  be  readily  incorporated 
if  desired.  These  polynomial  equations  are  given  in  Appendix 
B. 
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Once  the  values  of  Ttth,  Rc,  and  yc  are  known  then: 


and  the  experimental  value  of  chamber  stagnation  temperature 


T«-  is  determined  from 
cexp 


T*. 

rexp 


( 


fs. isJBJaV 

'*ave  +  "fV 


where  ^  *  throat  area 


CD  =  .97 

Finally,  the  efficiency  is  computed 


n 


AT 


Ttexp. 

Ttth  "  T 


Hard  copy  printouts  of  all  data  is  then  provided  by  the 
system  printer. 
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APPENDIX  A 


1 

2 

701 

705 

709 

721 

722 

723 

724 

"  :D700" 
" :D701" 

t«  ,  fji  •• 


LIST  OF  HP-IB  ADDRESSES 

CRT  (internal) 

Printer  (internal) 

2631B  Line  Printer 
7225  Plotter 

3497A  Data  Acquisition/Control  Unit 
HP-85  Controller 

3456A  High  Accuracy  Digital  Voltmeter 

6942A  Multiprogrammer 

3437A  High  Speed  Digital  Voltmeter 

HG-78  Scanivalve  Interface  #1 

HG-78  Scanivalve  Interface  #2 

Disc  Drive  #0 

Disc  Drive  #1 

Internal  Tape  Drive 
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APPENDIX  B 


PEPCODE  REGRESSION  EQUATIONS  FOR  HTPB  AND  PMM 

HTPB  =  Fuel  Type  fSTOICH  =  *073  PP  *  '0332  #:n/in3 

0  =  -z — -  Pp  =  40-60  psia 

“STOICH  U 

Tt  =  -933.3  +  48930  +  2O302  -  89903  -  205. 906  +  1.1925?. 

TH  A 

+  1.014  (10~4)TA2  -  1.84  (10'8)Ta3  -  1.138Ta  +  .3864TA0 

-  7.67(1O_5)Ta20 

=  1.40162  -  .277230  +  .148302  -  .O244803  +  .00523®6 

-  4 .994 (10~5)Ta  +  4 .255 (10~9)Ta2  -  2 . 36 ( 10 "1 3 ) TA3 

+  7.7585  (10'"5)Ti?  -  3.1138  (1O'5)T.  02  -  1 .885  ( 1 0  ”  9 )  T  2  D 

AAA 

Rc  =  52.966  +  3.5160  -  9.55t>2  +  5.7103  +  .  2739  06 

-  8.08 ( 1 0  -  4  ) T A  +  1 . 36 ( 10_7)TA2  -  5 . 99 ( 1 0" 1 1 ) ?A3 

+  1.0208  ( 10~3 ) T. 0  -  1.872  (10"4)T. J2  +  2 . 18 5 ( 1 0~7 ) T  2 0 


I 


n 


■» 


PMM  =  Fuel  Type 


f STOICH 
p  =  40-60 


0.12 


psia 


pp  =  .0426  #m/in' 


s»  1642  -  59020  +  14688  0^  -  "74780  +  4.05T^4i 

TH 

-  .18Ta  -  4 . 68Ta02  +  1.28TAi4  -  3.259  (10  5)TA24>2 

YC  =  1.34  -  .02650  -  9.126  (10"6)Ta  -  .1934>2  +  .  132903 

-  4.205  (10~5)fA*  +  7.429  (10~5)Ta?2  -  2.72  (10  5^A*4 

-9  —  22 

+  2.09(10  )Ta  t> 

Rc=  47, 1569+0  [27. 54^-51  2. 9108- p  (41. 242-o[l9- 545-11. 567^])'] 

+  (10'6)TA!l2403  -  3.2029(10"3)TA2; 
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APPENDIX  C 


PROGRAM  LISTINGS 


10  !  HELP  FOR  HEATER 
20  OPTION  BASE  i 

30  COM  ¥(100,4)  ,M9,N0«m  ,P0,V9,B9 
40  CRT  IS  2 

SO  DISP  "The  HEATER  program  contains  project  control  an 
d  “ 

60  DISP  "  DATA  ACQUISITION  and  DATA  REDUCTION  programs 
for  the  VITIATED  AIR  HEATER," 

70  DISP  “PRESS  '  SET  UP  *  KEY  TO  MAKE  PRELIMINARY  AIR  FLO 
U  CALCULATIONS." 

80  DISP  "WHEN  YOU  ASE  SATISFIED  WITH  YOUR  CALCULATIONS 
YOU  CAN  FLOW  THE  AIR  SYSTEM  TO  CHECK . “ 

90  DISP 

100  DISP  "PRESS  'FLOW'  KEY  TO  COLD  FLOW  THE  AIR  TO  MEAS 
URE  Pti  AND  Ttl  AND  COMPUTE  «:dot,AIR.“ 

110  DISP  “IF  YOU  ARE  NOT  SATISFIED  WITH  Ptl  SETTING  OR 
ACTUAL  AIR  FLOW  THEN  'FLOW'  IS  RUN  AGAIN." 

120  DISP  "AFTER  'FLOW'  IS  EXECUTED  THE  'START'  OPTION  K 
EY  IS  PRESSED  TO  START  HEATER  SEQUENCE." 

130  DISP  "'START'  BRINGS  HEATER  UP  TO  DESIRED  TEMPERATU 
RE  ." 

140  DISP  "WHEN  DESIRED  TEMP  IS  ATTAINED  THE  'HTRMNT'  HE 
ATER  MAINT  ROUTINE  IS  RUN  EVERY  20  SECS." 

ISO  DISP  "THE  HEATER  IGNITER  IS  DE-ENERGIZED  AFTER  8  SE 
CONDS  OF  OPERATION  OR  WHEN  700  DEG  IS  " 

160  DISP  "REACHED  ,  ALSO  IF  THE  AIR  TEMP  DROPS  BY  50  DE 
G  IN  5  SECONDS  DURING  START  UP  THE  IGN" 

170  DISP  "  IS  SHUT  OFF  AND  THE  HEATER  IS  SECURED." 

180  DISP  "WHEN  THE  AIR  TEMPERATURE  IS  WITHIN  150  DEGREE 
S  OF  THAT  DESIRED,  THE  HEATER  MAINT  " 

190  DISP  "  ROUTINE  IS  CALLED.  THIS  PROGRAM  EXECUTES  EV 
ERY  20  SECONDS  UNTIL  ANOTHER  OPTION" 

200  DISP  "  IS  SELECTED,  MEASURING  AND  RECORDING  THE  AIR 
TEMP  .  “ 

210  DISP  "  THE  SOLID  FUEL  RAMJET  OPTION  KEY  'SFRJ'  IS  P 
RESSED  TO  ENERGIZE  THE  SFRJ  WITH  OR  " 

220  DISP  "WITHOUT  THE  HEATER  IN  OPERATION.  THIS  OPTION 
RUNS  FOR  A  PRESET  TIME  BEFORE  SECURING" 

230  DISP  "  'SFRJ2'  THE  DATA  REDUCTION  PROGRAM  IS  AUTOMA 
TICALLY  CALLED  AND  EXECUTED  AFTER  'SFRJ'" 

240  CRT  IS  1 

2S0  CHAIN  "HEATER" 

260  END 
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4 

>1 


I 


10 

!  HEATER 

20 

OPTION  BASE  1 

30 

COM  V( 100 .4) ,M9 , 

NO*  C 1 1 

,P0  ,V9,B9 

40 

ON 

KEY* 

1 ,  "SET  UP"  GOSUB  150 

50 

ON 

KEY* 

2,  “FLOW" 

GOSUB 

790 

60 

ON 

KEY* 

3  > "SFR  J " 

GOTO 

2271 

70 

ON 

KEY* 

7, “TJTC" 

GOSUB 

1220 

80 

ON 

KEY* 

5 , "HELP  " 

GOSUB 

1240 

90 

ON 

KEY* 

6, "START 

"  GOSUB  1270 

100  ON  KEY*  4, “SHUTDOWN-  GOSUB  2140 
110  CLEAR  @  KEY  LABEL 
120  DISP  "SELECT  OPTION* 

130  GOTO  130 

140  !  Initialize  the  constants  as  needed  to  do  the  prel 
minar  y 

ISO  !  calculations  for  setting  up  the  air  heater. 

160  DISP  "WILL  YOU  BE  RUNNING  THE  HEATER  ONLY?  (Y/N>* 
170  INPUT  Y* 

180  IF  Y$="N"  THEN  GOTO  2880 
190  ! 

200  OUTPUT  709  j  "D04 , 4 , S , 6 , 7" 

210  LET  V9-280  !  STORAGE  TANK  VOLUME 
220  G-l.4  !  AIR  GAMMA 

230  DISP  “ENTER  STORAGE  TANK  GAGE  PRESSURE  (PSIG)" 

240  INPUT  B8 

250  LET  C= . 97  !  DISCHARGE  COEF . 

260  LET  Gi-32.2  !  GRAVITY 

270  LET  D=1 . 5  !  INCHES 

280  D1-.235  !  Dstar  choked  venturi 

290  LET  R=»53 . 3  !  GAS  CONSTANT  FOR  AIR 

300  CLEAR 

310  DISP  “ENTER  DESIRED  MASS  FLOW  RATE  (#/SEC>  M9" 

320  INPUT  M9 

330  PRINT  "  FOR  DESIRED  MASS  FLOW  RATE  OF  " jM9;“  */SEC . 
•• 

340  !  DISP  "ENTER  NOZZLE  DIAMETER  (Di-COLD  AIR=D*>" 
350  !  INPUT  Di 

360  DISP  "ENTER  COLD  AIR  TEMPERATURE  ( DEG=F )  Tl" 

370  INPUT  Tl 
380  LET  Ti=Tl+460 

390  DISP  "ENTER  DESIRED  HOT  AIR  TEMPERATURE  (DEG-R)  T8" 
400  INPUT  T8 

410  DISP  "ENTER  BAROMETRIC  PRESSURE  SETTING  (INCHES  HG > 
<29. 92-STANDARD  SEA  LEVEL)" 

420  INPUT  B9 
430  B9-.4913*E<V 

440  !  B9  WILL  BE  USED  TO  CONVERT  ABS  PRESSURES  TO  GAGE" 

450  !  NOW  COMPUTE  VISCOSITY  FOR  LATER  USE" 

460  LET  U= . 0000000229*T8*1 . 5/ < T8+ 198 . 6 > 
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470  ! 

480  !  COMPUTE  FUEL  AIR  RATIO  FROM  REPCODE  EXCERT 
490  LET  F ■ ( T 8-550  > /7800Q 

S00  !  This  F  was  for  air  at  approx.  70  degress  farenhei 
t . 

S10  !  compute  ethylene  flow  rate 
S20  LET  M8=F*M9  !  MOOT  FUEL 
530  M7»96/28*M8  !  OXYGEN  FLOW  RATE 
S40  V«183720*M7  !  OXYGEN  SET  VOLTAGE. 

550  V**VAL*  <  V ) 

56<’  OUTPUT  709  ;"AQ3,0,"4V1 
570  CLEAR 

580  DISP  USING  590  j  M8;826.06*MB 

590  IMAGE  "SET  ETHYLENE  FLOW  RATE*"  , ID .5D ,  "  #/SEC"  ,5X  ,?> 
D.3D,"  SCFM" 

600  DISP  USING  610  }  M7;722.S8*M7 

610  IMAGE  “SET  OXYGEN  FLOW  R ATE= “ , ID . 5D , "  #/SEC " , 5X , 2D . 
3D,"  SCFM" 

620  PRINT  "OXYGEN  VOLTAGE*  M;V;"  MILLIVOLTS" 

630  LET  A*PI/4*DiA2 

640  LET  F 1*SQR ( G* < 2/ < G+l ) ) A  < (G+l )/<G-i)>) 

650  PRINT  “F < GAMMA )=>"  jF  1 
660  LET  Q9»SQR(T1*R/G1) 

670  LET  P6*M9*Q9/(C*A*F1>  !  PSIA 

680  LET  P5*P6-B9  •  PSIG 

690  DISP  USING  700  j  P5 

700  IMAGE  "SET  PT  TO  “,6D,"  PSIG" 

710  LET  T5»<B8-P5*2>*V9*i44/<R*Tl*M9> 

720  DISP  USING  730  ;  TS 

730  IMAGE  "RUN  TIME  IS  ",5D.D,"  SECONDS’ 

740  COPY 
750  CLEAR 

760  DISP  "SELECT  OPTION" 

770  KEY  LABEL 

780  RETURN 

790  !  FLOW  OPTION 

800  OUTPUT  709  ; “DQ4 , 4 ,5 , 6 , 7"  !  AIR  VALVE  CLOSED 
810  CLEAR 
820  WAIT  2000 

830  !  THIS  SEGMENT  FLOWS  THE  SYSTEM  COLD  TO  DETERMINE  A 

CTUAL  MASS  FLOW  RATE 

840  OUTPUT  709  >"AC5"  !  MEASURE  PT1 

850  OUTPUT  722  ; "HSMO 02SW2Z0SO1L1 . 1STIS0F 1R1T3QX1 • 

860  GOSUB  2230 
870  V0*V 

880  PRINT  "VO*", VO 

890  OUTPUT  709  ;"DC4,6"  !  AIR  VALVE  OPENS 
900  WAIT  1000 

910  OUTPUT  722  ; "HSM002SW2Z0SO1L1 . 1STIS0F1R 1T3QX1 " 

920  GOSUB  2230 
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930  PRINT  "V=" ,V 

940  P»79.434*(V-VQ>  »  VOLTS  TO  PSIG 
950  PRINT  "PT1»",P 

960  OUTPUT  709  j"AC6"  !  MEASURE  TT1 

970  OUTPUT  722  } "HSW2SM0O2ZOSO1L1 . ISTISOFiRiTiQXi " 

980  GOSUB  2230 

990  T«41S7S*V+S36.6  !  VOLTS  TO  DEGREES  R 
1000  OUTPUT  709  ;"D04,6"  !  CLOSE  AIR  VALVE 
1010  !  CONVERSIONS  ON  T  AND  P  NEED  TO  BE  DONE. 

1020  Q9«SQR(R*T/G1) 

1030  M»C*A*F1*<P+B9)/Q9  «  ACTUAL  MA  FLOW 

1040  DISP  "ACTUAL  MASS  FLOW  RATEE  IS  ■*  ;Mj"  *n/s«c  .  M 

1050  M8*F*H  9  M7*96/28*M8 

1060  V=183720*M7  !  OXYGEN  SET  VOLTAGE 

1070  OUTPUT  709  > "A03 , 0 , "&VAL4 < V ) 

1080  P?»M9*Q9/<F1*C*A>  !  PSIA 

1090  P«P7-P-B9  »  PRESSURE  INCREMENT 
1100  PRINT  "  DELTA  PTi=",P 
1110  DISP  USING  1120  j  P,M9 

1120  IMAGE  "CHANGE  Pti  BY  “,40,"  PSIG  TO  ACHEIVE  A  Mdot 
<  AIR )  OF  “  ,  D . 5D , "  #«qss/s«c." 

1130  DISP  “ARE  YOU  SATISFIED  WITH  THE  MASS  FLOW  RATE  A 
ND  DOME  PRESSURE  SETTING?  <Y/N>" 

1140  INPUT  Y* 

1150  IF  Y$= "Y "  THEN  GOTO  1180 

1160  Ti*T  !  P6*P7!  SAVE  LAST  PTi  AND  TTi 

1170  GOTO  790 

1180  DISP  "SELECT  OPTION* 

1190  M9=M  !  ACTUAL38 DESIRED  MDOT 

1200  KEY  LABEL 

1210  RETURN 

1220  BEEP  200,100 

1230  RETURN 

1240  !  HELP 

1250  CHAIN  "HELP" 

1260  RETURN 

1270  !  HEATER  START  UP  ROUTINE 
1280  OUTPUT  709  j*DC4,6“  »  AIR  ON 
1  ?9fl  n  PAD 

1300  OUTPUT  709  j"AC5"  !  Pti 

1310  OUTPUT  722  > "HSW2SM002ZOS01L1 . 1STIS0F1R1T 1QX1 " 

1320  GOSUB  2230 

1330  P»79.434*(V-V0)  !  VOLTS  TO  PSIG 
1340  OUTPUT  709  j"AC6"  !  Ttl 

1350  OUTPUT  722  ; "HSW2SM002Z0SO1L1 . ISTISOFIRITIQXI " 

1360  GOSUB  2230 

1370  T*41575*V+S36.6  «  T  COLD 

1380  Q9-SQR<R*T/G1> 

1390  M«C*A*F1*<P+B9)/Q9  !  ACTUAL  MDOT  AIR 


1400  M8=F*M  9  M7=96/28*M8 

1410  V=18372Q*M7  !  02  SET  VOLTS 

1420  OUTPUT  709  ; "A03 , 0 , "6VAL* < V ) 

1430  OUTPUT  709  >"AC7“  !  T  hot 

1440  OUTPUT  709  >"DC4,7“  »  IGNITER  ON 

1450  OUTPUT  709  },‘DC4>4,5H  !  02  4C2H4  ON 

1460  ON  TIMER*  1,8000  GOSUB  1670 

1470  T3«T  »  T  LAST 

1480  WAIT  5000 

1490  OUTPUT  722  > "HSW2SM002Z0SQ1L1 .1STIS0F1R1T1QX1" 
1500  GOSUB  2230 

1510  T2*1002.3*V+457.53  !  Thot 

1520  T1-T2-T3+50 

1530  IF  T2>700  THEN  GOSUB  1720 

1540  IF  T1 <  0  THEN  1760 

1550  PRINT  HTshot=" ,T2 

1560  !  CHECK  TO  SEE  IF  HEATER  IS  WITHIN  DESIRED  LIMITS 

1S70  IF  T2+150-T8>0  THEN  1840 

1580  !  HEATER  IS  NOT  YET  NEAR  OPERATING  RANGE 

1S90  OUTPUT  709  >"AC7“ 

1600  T3=T2  !  LAST  T  HOT 
1610  WAIT  4000 

1620  OUTPUT  722  > "HSW2SMOQ2ZOS01L1 . 1STIS0F1R 1T1QX1 " 
1630  GOSUB  2230 

1640  T2=10Q2.3*V+4S7.S3  !  Thot 
1650  Tl=T2-T3+50 

1660  GOTO  1540  !  READ  T  HOT  AGAIN 
1670  !  SHUT  OFF  IGNITER 
1680  OUTPUT  709  > "DQ4 ,7" 

1690  PRINT  “IGNITER  SHUT  DOWN-5  SECS . “ 

1700  OFF  TIMER*  1 
1710  RETURN 

1720  OUTPUT  709  ;"DQ4,7“ 

1730  OFF  TIMER*  1 

1740  DISP  “  T : HOT  >  7QU  R  IGNITER  OFF " 

1750  RETURN 

1760  !  TEMP  GRAD  IS  NEG  TIME  FOR  SHUT  DOWN 
1770  OUTPUT  709  )"D04,7" 

1780  OUTPUT  709  >“004,4,5" 

1790  OUTPUT  709  >"A03,0,0" 

1800  OUTPUT  709  >“D04,6"  !  AIR  OFF 

1810  DISP  "  HEATER  SHUTDOWN-NEG  TEMP  GRADIENT* 

1820  OFF  TIMER*  1 
1830  STOP 

1840  !  HEATER  IS  NEAR  OPERATING  POINT 

1850  OFF  TIMER*  1 

1860  PRINT  "  HEY  I'M  CLOSE" 

1870  ON  TIMER*  3,10000  GOSUB  1940 
1880  DISP  "ENTER  HEATER  MAINT  ROUTINE’* 

1890  CLEAR 
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1900  DISP  “HEATER  RUNNING- 

1910  WAIT  500 

1920  KEY  LABEL 

1930  GOTO  1690 

1940  !  HEATER  MAINTENANCE 

1950  OUTPUT  709  ;"AC7" 

1960  OUTPUT  722  * "HSM002SW2ZQSO1L1FL0 . 0 1STIS0F1R 1T1QX1 • 

1970  GOSUB  2230  !  TAKE  READING 
I960  T2-1002. 3*0+457.53  ‘  Thot 
1990  PRINT  “Tihot iAIR»“)T2 
2000  T4-T8-T2  !  DELTA  Tjhot 
2010  IF  ABS(T4) >200  THEN  GOTO  2140 
2020  IF  ABS(T4 ) < 100  THEN  1940 
2030  IF  T4<  0  THEN  GOTO  2060 
2040  DISP  "  NOT  HOT  ENOUGH  “ 

2050  GOTO  2080 
2060  DISP  *  TOO  HOT  “ 

2070  GOTO  2080 
2080  F4*T4/780Q0 

2090  F=»F+F4  !  NEW  FUEL  AIR  RATIO 
2100  M8=F*M9  8  M7»96/28*M8 
2110  V-183720*M7 

2120  OUTPUT  709  ) "A03 , 0 , "4VAL*  < V ) 

2130  RETURN 

2140  !  HEATER  SHUTDOWN 
2150  OUTPUT  709  j“A03,0,0* 

2160  OUTPUT  709  ; “D04 , 7, 4 ,5“ 

2170  PRINT  "HEATER  SHUTDOWN" 

2180  PRINT  -  LAST  T:hot=>"*T2 
2190  WAIT  5000 

2200  OUTPUT  709  ; “D04 , 4 ,5 , 6 , 7 " 

2210  OFF  TIMER#  3 
2220  RETURN 

2230  !  SRQ  SERVICE  ROUTINE 
2240  STATUS  7,1  j  Z9 
2250  ENTER  722  ;  V 

2260  ENABLE  INTR  7  >8  8  CLEAR  722  8  RETURN 
2270  !  SOLID  FUEL  RAMJET  OPTION 

2280  OUTPUT  709  ; "D04 , 8 , 9 , 1 0 , 1 1 , 12 , 13“  !  AIR  DUMP 
2290  OFF  KEY#  4  8  ON  KEY*  4 ,  "SHUTDOWN  **  GOSUB  30S0 
2300  IF  N0*-"N-  THEN  2330 
2310  OUTPUT  709  i"DC4,9,13" 

2320  DIM  00*[901 

2330  O0*»“WF,3.i,0,3.2,0,3.3,0T,SF,2,3,i, . 001 , 12T,MI ,2 , 
1000T , WF , 3 , 999T , AC , 3T , WF , 4 . 2 , 1T“ 

2340  OUTPUT  723  *00* 

2350  OUTPUT  723  > "OP , 4 , . 03ST * 

2360  CLEAR  722  8  OUTPUT  722  * "HSM002SW2Z0SOiLiFL0 . 01STI 
S0F1R1T3Q" 


2370 

ON" 

2380 

2390 

2400 

2410 

2420 

2430 

2440 

24S0 

2460 

2470 


OUTPUT  709  ;  "DC4 , 14."  8  DISP  -READY  FOR  SFRJ  IGNITI 

OUTPUT  709  >"AC13"  '  Pcha«b*r 
OUTPUT  722  >"Xi"  8  GOSUB  2230 
U9»U  INITIAL  Pc  READING 
OUTPUT  709  )“AC14“  !  IGNITION 
OUTPUT  722  ;"Xi"  8  GOSUB  2230 
IF  ABS(V) <3  THEN  2420 
1  BEGIN  SFRJ  START  UP 
SETTIME  0,0 

OUTPUT  709  j "DC4 ,8,11,12" 

WAIT  10004T6 


2480  OUTPUT  709  j"D04, 12,11" 

2490  OUTPUT  709  >"AC13"  !  Pchanbtr 
2S00  OUTPUT  722  >"X1“  8  GOSUB  2230 
2510  IF  U> .75*<P0/1S1 .5)  THEN  2540 
2520  CLEAR 

2530  DISP  “NO  IGNITION  TRY  AGAIN"  8  GOTO  2420 
2540  !  BEGIN  MEASUREMENT  SEQ. 


2550  ON  TIMER*  1,T6*1000  GOSUB  3050 

2560  OUTPUT  722  > "HSMQ02SU2Z0SO1L1FLQ . 0 1STIS0F1R 1T3Q" 
2570  FOR  J=1  TO  100 

2580  OUTPUT  709  ;"ACS"  8  OUTPUT  722  j"Xi"  8  GOSUB  2230 
2590  P-79.434*<V-U0>  8  OUTPUT  709  }"AC6"  8  OUTPUT  722  } 
"Xi"  •  GOSUB  2230 

2600  T-415754V+536.6  8  U< J, 1)»C*A*F1*(P+B9>/SQR(R*T/G1) 


2610  OUTPUT  709  ;"AC12"  8  OUTPUT  722  >“Xi"  8  GOSUB  2231 


2620  T3-4 1575*0+536. 6  8  OUTPUT  709  j"AC7“  8  OUTPUT  722 
j-Xl"  8  GOSUB  2230 

2630  V<  J,4>=»1002.3*U+457.S  8  T2«V(J,4)  !  T:hot 
2640  OUTPUT  709  j "AC8"  8  OUTPUT  722  j"Xi"  8  GOSUB  2230 
2650  P«79.434*<U-U0>  8  OUTPUT  709  >"AC9"  8  OUTPUT  722  > 
"XI"  8  GOSUB  2230 

2660  T»41S7S*U+536.6  8  M3=C*A*F i* < P+B9) /SQR <R*T/Gi > 

2670  M4«U(J,1)~M3  8  V(J,3>*TIME 
2680  IF  N0*="N"  THEN  M4»0 

2690  OUTPUT  709  j"ACil"  8  OUTPUT  722  ;"Xi"  8  GOSUB  2230 


2700  P4«79.434*<V-V0>  8  OUTPUT  709  ;"AC10"  8  OUTPUT  722 


;"X1"  8  GOSUB  2230 

2710  T4*41575#0+536 . 6  8  V< J , 2 ) » <M3*T4+M4*T3> /V( J , 1 ) 
bar 


2720 

2730 

2740 

2750 

2760 


PRINT  "J»“ ;J, "TIME*"}M<J,3> 
PRINT  MTh«ad»",T4 
PRINT  "Ph*ad*",P4 
PRINT  "T :hot=",V<J,4> 

PRINT  :  dot  AIR»",V(J,i) 


»  T 
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2770  PRINT  : do t  PRI-“,M3 

2780  PRINT  "T  jbypass*’*  ,T3 
2790  PRINT  "T:BAR-",V( J,2> 

2800  NEXT  J 

2810  O0»SPOLL(723>  ‘  WAIT  FOR  HEM  INTR 
2820  IF  00064  THEN  2810 

2830  DISP  “ MUL T I P R OGR AMMER  HAS  P:chanb«r  data  “ 

2840  SEND  7  ;  UNL  MLA  TALK  23  SCG  12 
28S0  ENTER  723  ;  018  DISP  “  SLQT*=“>01 
2860  OUTPUT  723  ; "DC , 3T , WF , 4 . 2, OT , OP , 4, OT" 

2870  CHAIN  “SFRJ2" 

2880  DISP  “WILL  YOU  BE  RUNNING  THE  SOLID  FUEL  RAMJET?  < 
Y/N)  “ 

2890  INPUT  Y< 

2900  IF  Y*»"Y“  THEN  GOTO  .940 
2910  CLEAR 

2920  DISP  “TURBOJET  COMBUSTOR  OPTION  IS  NOT  AVAILABLE  A 
T  THIS  TIME.  SELECT  ANOTHER  OPTION" 

2930  GOTO  ISO 
2940  !  SFRJ  OPTIONS 
2950  CLEAR 

2960  DISP  “ENTER  DESIRED  CHAMBER  PRESSURE  (PSIA)“ 

2970  INPUT  PO 

2980  DISP  “ENTER  THE  DESIRED  BURN  TIME  IN  (SECONDS)- 
2990  INPUT  TS 

3000  DISP  “ENTER  THE  IGNITION  TIME  IN  (SECONDS)” 

3010  INPUT  T6 

3020  DISP  “WILL  THE  BYPASS  OPTION  BE  USED?  (Y/N)“ 

3030  INPUT  N0« 

3040  GOTO  190 

30SO  !  SFRJ  SHUTDOWN 

3060  OUTPUT  709  > “D04 , 8 ,9 , 1 1 , 12 , 13“ 

3070  OUTPUT  709  ;’’DC4,iOH  !  N2  PURGE  ON 
3080  WAIT  S000  8  OUTPUT  709  } "D04, 10" 

3090  GOSUB  2140  !  HEATER  SHUTDOWN 

3100  CLEAR  8  DISP  "SFRJ4  HEATER  SHUTDOWN  “  8  KEY  LABEL 
3110  RETURN 


10  !  SFRJ2  DATA  REDUCTION 
20  OPTION  BASE  1 

30  COM  V(1Q0,4>  >M9>N0fCU  ,P0,V9,B9 
40  SHORT  P<1000> 

SO  GCLEAR 
60  CLEAR 

70  Gi-32.174  !  qc 

80  C=* .  97  !  DISCHARGE  COEF  . 

90  DISP  “SFRJ  DATA  REDUCTIONS" 

100  DISP  “ENTER  INITIAL  GRAIN  WEIGHT* (GRAMS > " 

110  INPUT  U1 

120  DISP  “ENTER  FINAL  GRAIN  WEIGHT* ( GRAMS ) " 

130  INPUT  W2 

140  DISP  “ENTER  GRAIN  LENGTH=(CM> " 

ISO  INPUT  L2 
160  L2*L2/2 . 54 

170  DISP  "ENTER  INITIAL  GRAIN  D I AMETER* < CM ) " 

180  INPUT  DO 
190  D0*D0/2.S4 

200  DISP  “ENTER  FINAL  AFT  DIAMETER*(CM) " 

210  INPUT  D2 
220  02*02/2 . S4 

230  Al«PI*i .2*2/4  !  AREA  THROAT 

240  DISP  "ENTER  FUEL  TYPE*(HTPB/PMM> " 

2S0  INPUT  Fl$ 

260  IF  F1*=“HTPB"  OR  F1*»“PMM"  THEN  290 

270  DISP  “ONLY  FUEL  TYPES  AVAILABLE  ARE  <HTPB>  OR  <PMM> 

N 

280  GOTO  240 

290  W3=  < U1-U2 ) /4S3 . 6  !  DELTA  WEIGHT  IN  #n 

300  !  ENTER  Pc  DATA 

310  SEND  7  j  UNL  MLA  TALK  23  SCG  5 

320  FOR  1*1  TO  997  STEP  4 

330  ENTER  723  USING  "#,K,K,K,K“  ;  P ( I > , P ( I+l ) , P ( 1+2 ) , P ( 
I-t-3) 

340  NEXT  I 
3S0  ENTER  7 
360  • 

370  MAT  P*(S)*P 

380  FOR  1*1  TO  1000 

390  P ( I )*15 . 87*P ( I )♦ . 22  !  PSIG 

400  NEXT  I 

410  Z0*AMAX(P) 

420  JO-AMAXROW  8  T0«J0*.06 

430  PRINT  "MAX  Pc*MiZ0;“  AT  TIME*"jT0;“  SECONDS" 

440  FOR  I»J0  TO  1000 
450  IF  P  < I ) > . 2*Z0  THEN  490 
460  S2=I 

470  PRINT  "S2* " ; S2 
480  GOTO  500 
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490  NEXT  I 

S00  FOR  I=JO  TO  1  STEP  -i 
S10  IF  P<1)>.2*Z0  THEN  550 
520  Sl-I 

530  PRINT  "Sl="jSi 
540  GOTO  560 
550  NEXT  I 

560  S3-S2-S1  !  ELEMENTS  INCLUDED 
570  S4-.06*Si  B  S5=.06*S2 

580  S6-.06*S3  !  TIME  TBF 

590  REDIM  P (S3) 

600  S7*SUM<P> 

610  P1*S7/S3+B9  !  Pc  bar<PSIA> 

620  REDIM  P(1000) 

630  GCLEAR  8  LOCATE  20 ,120 ,20 >80  B  FRAME  B  SCALE  0,60,0 

,100 

640  LAXES  5,100,0,0,2,5,5 
650  FOR  1*2  TO  1000  STEP  2 
660  PLOT  I*.06,P<I> 

670  NEXT  I 

680  MOVE  S4 , 0  @  PLOT  S4,0  8  PLOT  S4,P(Si) 

690  MOVE  SS,0  t  PLOT  SS,0  9  PLOT  S5,P(S2) 

700  MOVE  JO*. 06,0  8  PLOT  JO*. 06,0  @  PLOT  J0*.06,P(J0) 
710  COPY 

720  !  FIND  M : DOT : AVE  AND  Tbar i AVE 

730  FOR  1*1  TO  100 

740  IF  U<  2,3X84  THEN  770 

750  V4-I 

760  GOTO  780 

770  NEXT  I 

780  FOR  I=V4  TO  100 

790  IF  V ( 1 , 3  >  <  S5  THEN  820 

800  V5-I 

810  GOTO  831 

820  NEXT  I 

830  V6*0  B  T-0  8  T3=0 
840  FOR  I-V4  TO  VS 
850  V6-V6+V(I,1)  !  M ; DOT  AVE 
860  T3*T3+V< 1,4)  »  T  i  HO  T  AIR 
870  T*T+V< 1,2)  !  T i BAR  AVE 
880  NEXT  I 

890  V8-V5-V4  !  TIME  INT 

900  V6-V6/V8  •  M : DOT  AVE 

910  T3-T3/V8  !  T  •.  HOT  !  AIR  :  AVERAGE 

920  T-T/V8  !  T:bar  AVE 

930  PRINT  "Pc  sbar*"  ; P 1 ;  *•  PSIA" 

940  PRINT  “Tbf *" j  S6  ; "  SECONDS" 

950  PRINT  "« : AVE*" j V6 
960  PRINT  "T ! BARave*" > T 
970  PRINT  "T:hot:air*";T3 
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980  PRINT 

990  V7=W3/S6  !  FUEL  FLOW  RATE 

1000  X1=V7/V6  •  fuel/air  ratio 

1010  PRINT  "f=" ;Xi 

1020  IF  F1*»"HTPB"  THEN  GOSUB  1130 

1030  IF  F1*="PMM"  THEN  GOSUB  1280 

1040  GOSUB  1420 

1050  GOSUB  1450 

1060  T2»<P1*A1*C*F2/<U6+V7> > A2  »  Tcexp 

1070  NO* <  T2-T ) / ( Tl-T  >  !  EFFICIENCY 

1080  PRINT  "Tc: THEORY-" *T1 

1090  PRINT  "Tc:EXP  *  ";T2 

1100  PRINT  "T : BAR : AVE*" ;T 

1110  PRINT  "THERHAL  EFFICIENY  =  "}N0 

1120  END 

1130  !  HTPB  SUBROUTINE 
1140  F9= . 073  !  FiSTOICH 
1150  D9-.0332  !  DENSITY  #«/inA3 
1160  X-X1/F9 

1170  T 1--933 . 3+4893*X+203*XA2-899*XA3-2Q5 . 9*XA6+1 . 1925* 
T  + . 0001014*TA2- . 0000000184*TA3-1 . 138*T*X 
1180  Ti=Tl+ . 3864*T*XA2- .0000767*TA2*X 

1190  G*1 .40162-. 27723*X+. 1483*XA2-. 02448*XA3+. 00523*XA6 

- . 00004994*T+ . 000000004255*TA2-2 .36E-13*TA3 

1200  G*G+. 000077585*T*X-. 00 0031 13B*T*XA2-. 00000 000 1085* 

TA2*X 

1210  R-52 . 966+3 . S16*X-9 . SS*XA2+5 . 71*XA3+ . 2739*XA6- .0008 
08*T  + . 000000 136 *TA2 -5 .99E— ll*TA3 

1220  R=R+ . 0 0 1 0208*T*X- .000 1872*T*X A2+ . 0  0000 02185*X*T A2 

1230  PRINT  "PHI* “  ; X 

1240  PRINT  "Tc : th= " ; Ti 

1250  PRINT  "Rc:th="  jR 

1260  PRINT  "GAMMAc*" ;G 

1270  RETURN 

1280  !  PHM  SUBROUTINE 

1290  F9*. 12  !  FiSTOICH 

1300  D9-.0426  «  DENSITY  #n/inA3 

1310  X-X1/F9 

1320  T i-1642-5902*X+14688*XA2-7478*XA  3+4 . 05*T*X- . 18*T-4 
. 68*T*X A2+i . 28*T*XA4- . 00003259*TA2*XA2 

1330  G=1 . 34- . 0265*X- . 193*XA2+ . 1329*XA3- . 000009126*T- . 00 
004205*T*X+ . 00007 429*T*XA2- . 0000272*T*XA4 
1340  G-G+  ,0000000020 9*TA2*XA2 

1350  R-47 . 1569+X* < 27 . 54+X*(2 . 91Q8-X*<4i . 242-X*( 19 . 54S-1 
1 .567*XA3) ) ) ) 

1360  R-R+. 00000 1*T*< 12403- . 0032029*TA2) 

1370  PRINT  "PHI* " j X 
1380  PRINT  "Tc : th* " j Tl 
1390  PRINT  "Rc:th="jR 
1400  PRINT  "GAMMAc* " ; G 
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1410  RETURN 
1420  !  F  <  GAMMA : C ) 

1430  F2*SQR<G*G1/R*(2/<G+1 ) ) * ( ( G+l ) / ( G-l ) ) ) 

1440  RETURN 

1450  !  Ds FINAL: AVERAGE 

1460  DB»SQR<4*U3/<PI*D9*L2>+D0A2)  !  D : FINAL  AUE  . 
1470  R8»<D8-D0>/<2*S6>  !  AUE  BURN  RATE 
1480  D8«2. 54*08  !  IN  CM. 

1490  R8»2.54*R8  !  IN  CM/SEC 
1S00  PRINT 

1510  PRINT  "d : final i AUE»“ ; D8 
1520  PRINT  "AUEiburn  RATE=MjR8 
1530  RETURN 
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10  !  HELP 

20  !  BERT  HANSEN 

30  ! 

40  !  THIS  PROGRAM  INITIATES  THE  LIGHT  DIFFRACTION  EXPER 
IMENT. 

50  !  IT  ALSO  PROVIDES  HELP  WITH  THE  DOCUMENTATION. 

60  • 

70  OPTION  BASE  i 
80  CLEAR 

90  COM  Qit [ 121 , Q2* (121 , Ti*  C  20 ] , Di$ [ 12 ] ,  D2* [ 12] 

100  DISP  “ENTER  DATE : TIME ( MMDDHHMMSS > “ 

110  INPUT  Ti* 

120  OUTPUT  709  j"TD“4T11 
130  OUTPUT  709  ;"TD" 

140  Q*= "DATA “  !  TYPICAL  NAME! 

150  CLEAR 

160  F=0  •  FLAG  RESET 
170  CRT  IS  1 

180  DISP  "ENSURE  THAT  THE  FOLLOWING  EQUIPMENT  IS  TURNED 
ON.  " 

190  DISP  "  1 . ALL  HP-IB  EQUIP." 

200  DISP  "  2 . VISICORDER " 

210  DISP  “  3 . RETICON  ARRAY." 

220  DISP  “  4. LASER" 

230  DISP  "  5. BATTERY  AND  IGNITION  CONTROL  CIRCUITS." 

240  DISP  "  6. LOW  PASS  FILTER" 

250  DISP  "The  following  programs  are  asssociated  with  t 
he  light  diffraction  experiment." 

260  DISP 

270  DISP  "  PRESS  CCONT3  WHEN  READY." 

280  IF  F=0  THEN  PAUSE 

290  DISP  "  l.ACQDTA-  acquires  data  from  both  diode  arra 
ys  and  stores  it  on  the  disc." 

300  DISP 

310  DISP  “  2 . CMBINE-c omb ines  raw  data  arrays  into  TWO  1 
024x2  arrays  and  datafiles." 

320  DISP 

330  DISP  "  3 . PLTDTA-p 1 o ts  the  raw  data  from  the  datafil 
e  as  VOLTAGE  VS.  DIODE  NUMBER." 

340  DISP 

350  DISP  "  PRESS  [CONTI  WHEN  READY." 

360  IF  F*0  THEN  PAUSE 

370  DISP  “  4 . REDUCE-per f orns  the  reduction  math  and  plo 
ts  NORMALIZED  INTENSITY  VS.  THETA." 

380  DISP 

390  DISP  "  5 . REDUC2-p 1 o t s  the  ideal  curve  from  disc.  Ca 
libration  data  can  also  be  plotted." 

400  DISP 

410  DISP  "  PRESS  [CONTI  WHEN  READY." 

420  IF  F-0  THEN  PAUSE 
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430  DISP  "  6,REDUC3--Conputes  D32  fro*  the  user  specific 
d  datafiieand  calcs.  Theta  bar  .  11 
440  DISP 

4S0  DISP  “  7.CALIB-  If  a.  run  is  to  be  a  calibration  cur 
ve  th  i  s  ■* 

460  DISP  "  proqraM  performs  a  curve  fit  to  get  a  snooth 
curve  for  conparison .  " 

470  DISP 

48C  DISP  “  8 . REDUC4-P lots  INTENSITY  US  THETA  BAR “ 

490  DISP  "using  D32  enterred  by  user  in  fornot  like  RED 
UC2  “ 

S00  IF  F=0  THEN  PAUSE 
510  IF  F-2  THEN  GOTO  580 

S20  DISP  "  DO  YOU  WANT  A  HARD  COPY  OF  THE  INSTRUCTIONS? 
< Y/N) " 

5 JO  INPUT  A* 

540  IF  A$=  "N "  THEN  GOTO  580 
S50  F*2  '  FLAG  SET 
560  CRT  IS  2 

570  GOTO  180  !  RETURNS  TO  PRINT 
580  CRT  IS  1 
590  CLEAR 

600  CHAIN  “ ACQDT A “ 

610  END 
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AUTOMATIC  CONTROL  AND  UA1A  ACOUlSHlUN  b'lbltM  hUK 
COMBUSTION  LABORATORY  APPL ICAT IONS t U )  NAVAL 
POSTGRADUATE  SCHOOL  MONTEREY  CA  B  d  HANSEN  OCT  82 

F/G  21/2 


UNCLASSIFIED 


100  •  ACQDTA 
110  !  BERT  HANSEN 
120  OPTION  BASE  1 

130  OUTPUT  709  "TOSOO"  !  TIME  SIGNAL  TO  VISICORDER 
140  COM  Q1 01121 , Q24I12) ,T14(201,D14(12] ,D24( 121 
150  SHORT  A<1032) ,8(1032) 

160  OUTPUT  709  }"TD" 

170  ENTER  709  }  T16 

180  PRINT  “DATE? TIME  IS  "jTil 

190  OUTPUT  709  "AR ,D04, 0 ,2, 3" 

200  OUTPUT  723  j "CC,2 ,3, 12 , 13T , WF ,3 . 1 , 0 ,3 .2, 0 ,3 .3 , QT, WF 
, 13. 1,0, 13. 2, 0,13. 3, OT" 

210  CLEAR 
220  CLEAR  722 

230  DISP  "ENTER  THE  FILENAME  OF  THE  DATA  FILES  TO  BE  CR 
EATED. (RAU1 ,RAW2)  (MAX  12  CHAR.)" 

240  DISP  "DATA  SHOULD  BE  STORED  ON  DRIVE  01  BY  TYPING 
( : D701 )  AFTER  FILENAME." 

250  DISP  "  BE  SURE  THE  FILES  NAMED  DO  NOT  ALREADY  EXIST 
ON  DISC" 

260  INPUT  D14,D26 
270  CLEAR 
280  M5-1024 

290  CREATE  014,1024,16 
300  CREATE  D24,1024,16 

310  PRINT  "  DATA  MILL  BE  STORED  ON  DISC  WITH  FILE  NAME 
AS  FOLLOWS . " 

320  PRINT  "NO  PARTICLES-" ,D14 
330  PRINT  "PARTICLES**"  D24 

340  DISP  "ENTER  THE  THRESHOLD  PRESSURE  TO  TRIGGER  DEVIC 
ES  (PSI)" 

350  INPUT  S8  !  DESIRED  PRESSURE 

360  DISP  "ENTER  TIME  DELAY  IN  SECONDS  FROM  THRESHOLD  PR 
ESSURES  (SECS.)" 

370  INPUT  T8 
380  TB«(T8-. 25)41000 
390  VO-38/151.5 
400  I 
410  CLEAR 

420  DISP  "  LA8ER  ON  , IGNITION  OFF" 

430  D34-D14 

440  GOSUB  840  t  CALL  MULTIPROGRAHMER 
450  G08UB  1180  !  STORE  DATA 

460  OUTPUT  709  "ACO"  !  CONNECT  3456  TO  IGNITION  FIRING 
LINE 

470  LOCAL  LOCKOUT  7 
480  R7-225 

490  OUTPUT  722  "HSM0028U2SO1L1S0F4R1T3QX1" 

500  GOSUB  1260  !  ENTER  READING 
510  R8>V 


S20  IF  R8<R7  THEN  GOTO  570 

em  I  pnklTTklt  ITTV  puprv 

540  DISP  “CONTINUITY  CHECKS  BAD  RECHECK  BEFORE  PROCEEDI 
NG,  THEN  PRESS C CONTI “ 

5S0  PAUSE 
560  GOTO  490 
570  CLEAR 

580  DISP  “CONTINUITY  IS  OK“ 

590  OUTPUT  709  j“DC4,2“  i  NITROGEN  ON 
600  DISP 

610  OUTPUT  709  j“AC2“  !  CONNECT  DVM  TO  FIRE  SWITCH 
620  DISP  *  ENSURE  UISICORDER  IS  SET  UP  TO  RUN  ON  PROPER 
SCALE  AND  LAMP  IS  ON* 

630  OUTPUT  722  “HSH002SU2Z0SO1L1FL0 . 0 1STIS0F1R 1T3Q* 

640  DISP  -STANDING  BY  FOR  IGNITION- 
650  BEEP 

660  OUTPUT  722  “XI*  •  GOSUB  1260  !  ENTER  DATA 

670  R9«ABS<U)  •  IF  R9<10  THEN  GOTO  660 

680  OO-TIME  «  OUTPUT  709  “004,0*  8  ON  TIHERO  1,5000  COS 

UB  1250 

690  OUTPUT  709  “ACi*  •  CLEAR 

700  OUTPUT  722  j“Xi“  R  GOSUB  1260  !  ENTER  DATA 

710  R9*ABS<V)  •  IF  R9<U0  THEN  GOTO  700 

720  WAIT  T8  •  01 -TIME 

730  GOSUB  840  I  1ULT1PROGRAMMER 

740  WAIT  3000 

7S0  OUTPUT  709  }“DO4,0“ 

760  LOCAL  7 
770  D34-D2* 

780  GOSUB  1180  !  STORE  DATA 
790  PRINT  “TIME  FIRE  TO  HULTI  CALL»“;Oi-QQ 
800  PRINT  “TIME  FIRE  TO  HEM  INTR»“j02-00 
810  PRINT  “DATA  STORED  ON  DISK  WITH  FILENAMES  <"jDl*j“> 
AND  < “ } D2f } “ ) “ 

820  CHAIN  “CMBINE" 

830  END 

840  1  MULTIPROGRAMMER  ROUTINE 

850  OUTPUT  723  j “SF ,2, 3, 1 , . 001 , 12T,SF, 12,3, 1 , . 101 , 12T* 
860  OUTPUT  723  j “WF ,3 . 1 , OT , WF,3 . 2 , OT, WF ,3 . 3 , OT, WF , 13 . 1 , 
OT ,WF, 13 . 2, OT , WF , 13 , 3 , OT* 

870  OUTPUT  723  } “WF,3,1031T,WF, 13,1031T,MI,2,1032T“ 

880  OUTPUT  709  >“DC4,3“  I  MARKS  THE  UISICORDER  TRACE 

890  OUTPUT  723  } “AC ,3T , WF , 4 . 2 , IT ,OP ,4 , 16. ST- 
900  K-8P0LL(723>  !  WAIT  FOR  MEM  INTERUPT 
910  IF  K<>64  THEN  GOTO  900 
920  02»TIME 

930  OUTPUT  709  *“D04,3" 

940  DISP  *  DATA  TAKEN  AND  STORED* 


950  SEND  7  ;  UNL  MLA  TALK  23  SCG  12 
960  ENTER  7  )  K1 

970  DZSP  "MEMORY  BOARD  HAS  DATA" 

980  DISP  "  SLOT*"" ;K1 

990  OUTPUT  723  j "DC ,3T ,UF ,4 . 2 , OT" 

1000  ENABLE  1NTR  7}8 

1010  SEND  7  j  UNL  MLA  TALK  23  SCG  S 

1020  FOR  1-1  TO  1025  STEP  8 

1030  ENTER  723  USING  "*  ,K  ,K  ,K  ,K  ,K  ,K  ,K  ,K  "  >  A<I),A<H-i>, 
A< 1+2) ,A<I+3),A< 1+4) ,A<I+S> ,A(I*6) ,A< 1+7) 

1040  NEXT  I 
1050  ENTER  7 

1060  DISP  "  ARRAY  ONE  DATA  ENTERRED" 

1070  OUTPUT  723  ) "MI , 12, 1032T" 

1080  SEND  7  j  UNL  MLA  TALK  23  SCG  5 
1090  FOR  1*1  TO  1025  STEP  8 

1100  ENTER  723  USING  "#,K ,K ,K  ,K  ,K  ,K  ,K  ,K "  j  B<I),B<I+1), 
B(I+2) ,B( 1+3) ,B< 1+4) ,B<I+5),B<I+6>,B<I+7> 

1110  NEXT  I 
1120  ENTER  7 

1130  DISP  "  ARRAY  TWO  DATA  ENTERRED" 

1140  MAT  A— A 
1150  MAT  B*-B 

1160  OUTPUT  723  } "CC,2,3, 12, 13T" 

1170  RETURN 

1180  !  THIS  ROUTINE  STORES  THE  RAW  DATA  ON  THE  DISK 
1190  ASSIGN*  1  TO  D3* 

1200  FOR  1*7  TO  1027  STEP  4 

1210  PRINT*  1  ;  A<I),A(I+1>,A(I+2),A<I+3),B<I),B<I+1),B 
(I+2),B(I+3) 

1220  NEXT  I 

1230  ASSIGN*  1  TO  * 

1240  RETURN 

1250  OFF  TIMER*  1  •  OUTPUT  709  j"DO4,0,2*  8  RETURN 

126*  I  8RQ  SERVICE 

1270  STATUS  7,1  •,  AO 

1280  ENTER  722  )  V 

1290  ENABLE  INTR  7)8  •  RETURN 


10  •  CMBINE 

20  f  COMBINES  RAW  DATA  FILES  INTO  A  SINGLE  DATA  FILE 
30  OPTION  BASE  1 

40  COH  Q10(12],Q2«(t2],T10C20],Dl«(12],D20(12I 

SO  D90*" :D70i " 

60  CLEAR 

70  ! 

80  SHORT  Y (1024,4) 

90  D30-D1* 

100  FOR  M-i  TO  2 
110  ASSIGN*  2  TO  D3* 

120  FOR  1*1  TO  1021  STEP  4 

130  READ*  2  i  Y<I>H)>Y(I>l,H)>Y<I+2,H>>Y<I+3,N),Y(I,N+2 
>,Y<I+i,M+2>,Y<H-2>M+2>,Y<I+3>M+2> 

140  NEXT  I 

ISO  ASSIGN*  2  TO  4 

160  D3t*D20 

170  NEXT  M 

ISO  Q34-D14 

190  FOR  N-1  TO  3  STEP  2 
200  ASSIGN*  1  TO  Q30 
210  FOR  N-1  TO  1017  STEP  8 

220  A*Y<N,M)  •  B-Y<N+1,M>  8  C*Y<N+2>N)  8  D-Y(N+3,H> 

230  E-Y<N+4,H>  •  F-Y<N+S,H>  «  G-Y<N+6>H)  8  H-Y<N+7,H> 
240  R-Y(N>H+1)  8  3*Y<N«-1  ,M*1>  8  T-Y<N*2,H*i>  8  U-Y(N+3> 
N+l)  8  g-Y(N+4,M+l) 

2S0  M-Y<N+5,H+1)  8  X*Y<N-«-6,H-H  >  8  Z“Y<N+7 ,H+1  > 

260  PRINT*  1  i  A,B,C,D,E,F,G,H,R,S,T,U,V>U>X,r. 

270  NEXT  N 

280  ASSIGN*  1  TO  * 

290  Q30-D2* 

300  NEXT  M 

310  D1SP  "  DATA  STORED  ON  DISK  ON  FILE  <">D10j">  AND  ( - 
i  D2*  j "  )  " 

320  PRINT 

330  PRINT  "DATA  FILES  HAVE  BEEN  REARRANGED  A  FOLLOWS:" 
340  PRINT  "PHOTO  ARRAY  *1  IS  ";D1« 

3S0  PRINT  "PHOTO  ARRAY  02  IS  "}D2* 

360  CHAIN  "PLTDTA" 

370  END 


10  !  PLTDTA  REVISION  1 
20  ! 

30  !  PLOTTING  RAW  DATA  FROM  ARRAY 

40  I 

SO  !  BERT  HANSEN 
60  CLEAR 

70  OPTION  BASE  1 

80  COM  Qlt(121,Q2t(121,Tlt(201,Dlt(12],D2t(121 
90  OUTPUT  709  j "TD" 

100  ENTER  709  j  Tit 

110  DISP  “  ENSURE  PLOTTER  IS  READY  FOR  PLOTTING  WITH  CL 
EAN  PAPER  AND  DESIRED  PEN, • 

120  DISP  "PRESS  (CONTI  WHEN  READY!" 

130  BEEP 
140  PAUSE 

ISO  DISP  "  ENTER  THE  FILE  NAME  OF  THE  DATA  TO  BE  PLOTTE 
D .  ( " jDlt j "  OR  "*D2t*"  )" 

160  DISP  "DRIVE  MUST  BE  SPECIFIED  AS  (NAMEiD701) " 

170  INPUT  Q3t 

180  DISP  "ENTER  DESIRED  PLOTTING  INTERVAL  OF  DIODE  ARRA 
Y  <2,4,6)" 

190  INPUT  J 
200  CLEAR 

210  DISP  "READING  DATA  FROM  DATA  FILE  ON  DISC  DRIVE 
■ 

220  PLOTTER  IS  70S 

230  SHORT  Y(1Q24,2) ,X< 1024) 

240  HS-1024 

2S0  REDIM  X<HS) ,Y<MS,2) 

260  Pit- "VOLTAGE" 

270  P2t- "DIODE  NUMBER" 

280  P3t— "VOLT AGE  VS.  DIODE" 

290  P4t- "PLTDTA  RESULTS" 

300  ! 

310  ASSIGN#  1  TO  Q3t 
320  FOR  1-1  TO  MS-7  STEP  8 

330  READ#  1  }  A,B,C,D,E,F,G,H,R,8,T,U,V,W,X,Z 

340  Y ( I , 1 ) -A  •  Y(I+1 , 1 )-B  •  Y<I+2,1>-C  8  Y(H-3,1)-D  B  Y 

<1+4, 1)-E  •  Y<I>S,i)-F 

350  Y<I+6,1)-G  t  Y< 1+7 , 1 >-M 

360  Y<I,2)-R  •  Y < I+l ,2)-S  •  Y(I+2,2)-T  9  Y<H-3,2>»U  •  Y 

< I+4,2)-V  •  Y<IfS,2>»W 

370  Y<H>6,2)-X  •  Y< 1+7 ,2>-Z 

380  NEXT  I 

390  M7-MAXAB<Y) 

400  DI8P  "DATA  INPUT  COMPLETE" 

410  LIMIT  0,250,0,180 
420  FRAME 

430  LOCATE  20,120,20,80 
440  FRAME 


100 


450  SCALE  0 ,M5,- . 1 ,M7+ . 1 
460  C3IZE  3,. 5,0 
470  DEC 
480  LDIR  • 

490  FXD  0,2 

500  LAXES  MS/20, .02,0,0,2,5,5 
510  HOVE  .65*M5, .93«<H7+. i) 

520  LABEL  P4t 

S30  MOVE  . 654MS , . 84*<M7+ . 1 > 

540  LABEL  P3t 
550  LABEL  Tit 
560  MOVE  MS/2 , - . 05 
570  LABEL  P2t 
580  LDIR  90 

590  MOVE  -< .1*M5) , .5*M7 
600  LABEL  Pit 
610  LDIR  0 
620  M«i 

630  GOSUB  760  !  PLOT  FIRST  ARRAY 
640  DISP  “  FIRST  ARRAY  PLOTTED  I* 

650  M«2 

660  GOSUB  760  !  PLOT  SECOND  ARRA492  HOVE  . 65*H5, . 74*(M7 

+  .1) 

670  DISP  "  SECOND  ARRAY  PLOTTED!" 

680  DISP  "  PLOT  COMPLETE!  * 

690  DISP  "DO  YOU  WISH  TO  PLOT  THE  OTHER  PHOTO  ARRAY  FIL 
E?  < Y/N) " 

700  INPUT  Nt 

710  IF  Nt-"Y"  THEN  GOTO  110 
720  CLEAR 

730  DISP  "LOAD  AND  RUN  < REDUCE)  TO  PROCEED  WITH  DATA  RE 
DUCTI0N8." 

750  END 

760  !  PLOT  DATA 
770  X<1)»1 

780  MOVE  X(1),Y(1,1) 

790  FOR  1-1  TO  MS  8TEP  J 
800  X(I)«I 

810  PLOT  X(I),Y(I,M> 

820  NEXT  I 
830  PENUP 
840  RETURN 
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10  !  REDUCE 
20  CLEAR  7 

30  D1SP  "REDUCE  LOADED" 

40  \ 

50  I  BERT  HANSEN 

60  i  THIS  PROGRAM  REDUCES  THE  ARRAY  DATA 
70  OPTION  BASE  1 
80  DIM  A(2,2),B(2),Y1(2) 

90  COM  Q30C12] ,0201 121 ,T10!201,D10tl2) ,D20C 121 ,0 
100  OUTPUT  709  >"TD" 

110  ENTER  709  }  T10 
120  CLEAR 

130  DISP  "DATA  REDUCTION  BEGINS" 

140  DISP  "  ENTER  THE  INDEX  OF  REFRACTION-? (AIR  N-i . 0 ,  U 
ATER  M-1.35)" 

ISO  INPUT  M 

160  DISP  "INSTALL  PAPER  AND  PEN  ON  PLOTTER  TO  PLOT  I<> 
VS  " 

170  DISP  "PRESS! CONTI  WHEN  DONE" 

180  PAUSE 

190  DISP  "ENTER  THE  FILE  NAME  OF  DATA  TO  BE  REDUCED." 
200  DISP  "INCLUDE  THE  DISC  LOCATION  ON  FILENAME  (  NAME: 
D701 ) " 

210  INPUT  AO 
220  D90-" : D701 " 

230  L-. 000000488  !  WAVELENGTH  IN  METERS 
240  SHORT  TH1024), 11(1024), 12(1024) 

2S0  MS-1024  !  CURRENT  NUMBER  OF  ACTUAL  READINGS 

260  D8-1.536S  I  MMETERS 

270  F8-S26  «  MMETERS 

280  D1-.02S  I  MMETERS 

28S  GCLEAR 

290  PLOTTER  IS  705 

300  LIMIT  0,250,0 , 180 

310  LOCATE  20,120,20,80 

320  FRAME 

330  SCALE  0,. 04, -.1,1 
340  CSIZE  3,. 5,0 
350  DEG 
360  LDIR  I 
370  FXD  3,1 

380  LAKES  . 002, . 1 , 0 , 0 ,2,2,5 

390  DISP  "DO  YOU  WANT  TO  PLOT  A  COMPARISON  CURVE  (  UNIC 
UR)?  (Y/N)" 

400  INPUT  NO 

410  IF  NO-"Y"  THEN  GOSUB  1140 

420  ASSIGN#  1  TO  AO 

430  FOR  I-t  TO  MS-7  STEP  8 

440  READ#  1  i  A,B,C,D,E,F,C,H,R,S,T,U,V,U,X,Z 

450  I2( I )-A  «  I2( I+l )-B  8  I2(I+2)-C  •  I2(I+3)»D  •  12(1+ 
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4>-E  •  X2U+S)-F  8  I2<I+6)»G  8  I2(I+7)*H 

460  I1(I)-R  8  I1<I+1)-S  8  I1(I*2)-T  8  Ii(I+3)«U  8  11(1+ 

4>-V  8  Ii(I+5)*U  8  II < 1+6  >«X  8  I1<I+7)«Z 

470  NEXT  X 

480  MAT  11*11-12 

490  Z9-AMAX(I1) 

S00  ZO-AMAXRQM 

510  PRINT  -MAXIMUM  OUTPUT  AT  DIODEO--JZO 

520  DXSP  -ENTER  THE  INCREMENT  AFTER  PEAK  < J-4START ,K«# 

STOP  FIT)- 

530  INPUT  J,K 

540  N0*K-J 

5S0  X1*0  8  X2*0  8  X3*Q  8  Yl-0  8  Y2*0 
560  FOR  N-J  TO  K 
570  T1<N)*N+61.46 
580  X1«X1+T1<N> 

590  X2-X2-»*T1<N>A2 
600  X3*X3^T1(N)*I1(N) 

610  Yt-Yl+IKN) 

620  Y2»Y2*I1(N)A2 
630  NEXT  N 
640  A(1,1)-NI 
650  A< 1 ,2)*X1 
660  A<2,1)*X1 
670  A(2,2)*X2 
680  Y1<1)»Y1 
690  Yl(2)-X3 
700  MAT  B«SYS<A,Y1> 

710  PRINT  "DELTA*" >B<1) 

720  PRINT  -SLOPE-" }B(2> 

730  A3-B< 1 ) 

740  R2-(N0*X3-X1*Y1>A2/(<N0*X2-X1A2>*<N0*Y2-Y1A2))  !  CO 

RRELATION 

750  PRINT 

760  PRINT  -CORRELATION  RA2**>R2 
770  MAT  U-U/A3)*I1 

780  !  II  ARRAY  IS  THE  FINAL  PLOTTING  NORMALIZED  ARRAY 
790  !  MOM  PLOT  II  VS.  T1 
800  FOR  1*1  TO  M5 

810  T1(I>-<I*.02S+D8)/<FB*M)  !  THETA 
820  PLOT  T1(X>,X1(I> 

830  NEXT  X 

840  PENUP 

850  MOVE  0,0 

860  MOVE  .02,-. 3 

870  LABEL  "THETA  (RADIANS)" 

880  MOVE  -.004,. 2  8  DEG  8  LDIR  90 
890  LABEL  -NORMALIZED  INTENSITY" 

900  LDIR  0  8  MOVE  .015,. 8 
910  LABEL  -INTENSITY  VS.  THETA" 
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920  HOVE  .015 ,.9 

930  LABEL  “REDUCE  RESULTS'* 

940  HOVE  .015 ,.7 
950  LABEL  Tit 
960  Q5t*"REDUC3" 

970  DXSP  “  IS  THIS  A  CALIBRATION  RUN?  (Y/N>“ 

Oflfl  TIMMIT 

990  IF  Nt*“N"  THEN  GOTO  1010 
1000  Q5t*"CALIB“ 

1010  DXSP  “ENTER  THE  NAHE  OF  THE  DATA  FILE  TO  STORE  THE 
DATA  < INTENSITY  VS. THETA )“ 

1020  DISP  “INCLUDE  THE  DISC  LOCATION  ON  FILENAHE  <  NAHE 
: D701 ) “ 

1030  INPUT  Q3t 

1040  CREATE  030,1024,16 

1050  PRINT  “INTENSITY  VS. THETA  STORED  IN  FILE-"iQ3f 

1060  ASSIGN#  2  TO  Q30 

1070  FOR  1*1  TO  HS 

1080  PRINT#  2  j  Tl( X) , Xl( I ) 

1090  NEXT  I 

1100  ASSIGN#  2  TO  4 

1110  DISP  “DATA  STORED  ON  DISK  IN  FILE  —  “jQ3t 
1120  CHAIN  QSt 
1130  END 

1140  !  COHPARE  WITH  UNICUR 
1150  A88IGN#  2  TO  "UNICUR “ 

1160  FOR  1*1  TO  100  STEP  2 

1170  READ#  2  }  Ti(I),Ii(I),Ti(I+i),Il<I+i> 

1180  NEXT  I 

1190  DISP  “ENTER  EXPECTED  VALUE  OF  Dbar  (23. 98E-6)ETC. “ 

1200  INPUT  D 

1210  FOR  1-1  TO  99 

1220  T1<I>-T1(I)*L/<PI*0> 

1230  PLOT  T1(I> ,11(1) 

1240  NEXT  I 
1250  PENUP 
1260  RETURN 
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10  »  CALIB 
20  OPTION  BASE  i 

30  CON  Q39(121 >020(121 >110 1201 >010(121 ,020112) ,M 
40  D99-"»D701" 

SO  •  CALIBRATION  PROCEDURE 
60  •  BERT  HANSEN 
70  CLEAR 
80  CCLEAR 

90  DISP  "ARE  YOU  RUNNING  CALIB  SEPARATELY?  <Y/N>" 

100  INPUT  NO 

110  IF  NO-"N"  THEN  160 

120  DISP  "ENTER  DATAFILE  NAME  AND  DRIVE  ADDRESS" 

130  INPUT  939 

140  DISP  "ENTER  INDEX  OF  REFRACTION  (AIR-1.0  ,UATER-1.3 
S>" 

ISO  INPUT  H 

160  DISP  "ENTER  THE  ADDRESS  OF  PLOTTER.  (CRT-1,  722SB-7 
OS)" 

170  INPUT  P 
180  PLOTTER  IS  P 

190  IF  P-70S  THEN  LIMIT  0,250,0,180 
200  LOCATE  20,120,20,80 
210  FRAME 

220  SCALE  0,10,-.1,1 
230  LAXES  .5, .2, 0,0, 2, S,S 

240  • 

2SQ  DIM  A(U,11),B(11),W(20>,Y1(U),Y2(11> 

260  SHORT  X( 1024) ,Y (1024) 

270  MAT  A-ZER 
280  MAT  B-ZER 
290  MAT  Yl-ZER 
300  MAT  U-ZER 
310  MS-1024 

320  DISP  "  INDEX  OF  REFRACTION  IS  -";M 
330  DISP  "USING  (  "jQSOj"  >" 

340  DISP  "A  POLYNOMIAL  REGRESSION  ROUTINE  OF  (J)TH  ORDE 
R  WILL  BE  PERFORMED.  " 

350  DISP  "ENTER  ORDER  OF  POLYNOMIAL (  REASONABLE).* 

360  INPUT  Ji 

370  A88IGM9  1  TO  Q39 

380  FOR  I-l  TO  MS  STEP  2 

390  READ#  t  i  X(1),Y(I) ,X( X+l ) , Y( I+l ) 

400  NEXT  I 

410  DI8P  "ENTER  Dbar  TO  USE  IN  DATA  REDUCTION  («.*.  23 
. 98E-6) " 

420  INPUT  D3 

430  L-. 000000480 

440  F-S26  I  FOCAL  LENGTH 

4S0  X6-D39PI/L 

46*  OR  I-l  TO  MS 
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470  X<I>*X<I)*K6 
480  PLOT  X(I) >Y(I) 

490  NEXT  I 

500  •  I  VS. THETA  BAR 
510  PENUP 

520  DISP  "ENTER  THE  MINIHUH  THETA  BAR  OF  INTEREST." 

530  INPUT  J 

540  J»M*J/<236*D3>-61.46 
550  K-l . 5365 

560  OISP  "ENTER  THE  MAXIMUM  THETA  BAR  OF  INTEREST  . <ba 
r«3.S>" 

570  INPUT  Ml 

S80  Mi«<Ml*L*F*M/<PI*D3)-K  >/ . 02533 

590  PRINT  "MIN  DIODE 

600  PRINT  "MAX  DIODE  «"}M1 

610  J2-J1+1 

620  FOR  1*1  TO  2*J1 

630  I 

640  S0»Q 

650  FOR  M*J  TO  Ml 
660  ! 

670  S0»X<M)*I+S0 
680  NEXT  M 
690  U<I)*S0 
700  NEXT  I 
710  ! 

720  FOR  1*1  TO  J2 
730  S0*0 

740  FOR  N-J  TO  Mi 

750  SO«SO+Y<N>*X<N)MI-l> 

760  NEXT  N 
770  Y1<I)-SQ 
780  NEXT  I 

790  NO-Mi-J  !  0  OF  DATA  PAIRS 
800  ! 

810  REDIM  A<J2>J2),B<J2),H<2*J1),Y1(J2),Y2(J2) 

820  A  < 1 , 1 ) *N0 

830  FOR  1*2  TO  J2  !  FIRST  ROW 
840  A<l,I)*W<I-i> 

850  NEXT  I 
86Q  « 

870  FOR  1*2  TO  J2  \  FIRST  COLUMN 
880  A< I , 1 >*W< 1-1 ) 

890  NEXT  I 
900  \ 

910  M*0 

920  FOR  1-2  TO  J2  !  CORE  10*10 
930  FOR  K-2  TO  J2 
940  A<I,K>-W<K+M) 

950  NEXT  K 
960  M'  1+1 


970  NEXT  I 

980  GOSUB  1300  !  CORRELATION 
990  • 

1000  HAT  B=*SYS(A,Yi ) 

1010  PRINT 

1020  PRINT  “POLYNOMIAL  COEFFICIENTS® : " 

1030  PRINT 
1040  MAT  PRINT  B 
1050  ! 

1060  !  COMPUTE  THE  CALIBRATION  CURVE 
1070  S«.l 

1080  FOR  1*1  TO  100 
1090  X<I)»I*S 
1100  80-0 

1110  FOR  N-l  TO  J2 

1120  S0-S0*8<N)*X<I>*<N-1> 

1130  NEXT  N 

1140  Y< I )*S0 

1150  PLOT  X(I),Y(I) 

1160  NEXT  I 

1170  DISP  "  COMPLETE  !" 

1180  DISP  “ENTER  THE  DESIRED  FILE  NAME  TO  STORE  DATA." 
1190  INPUT  Q3* 

1200  CREATE  034,100,16 
1210  ASSIGN*  2  TO  Q3* 

1220  !  STORE  SMOOTH  CAL  CURVE  ON  DISC. 

1230  FOR  1-1  TO  100  STEP  2 

1240  PRINT*  2  }  X(I),Y<I) ,X<I+1) ,Y(I+i) 

1250  NEXT  I 

1260  DISP  "DATA  STORED  ON  DISC" 

1270  ASSIGN*  2  TO  * 

1280  CHAIN  “REDUC2" 

1290  END 

1300  !  COMPUTE  CORRELATION  FACTOR 
1310  DI3P  "  CORRELATION  " 

1320  SHORT  X4(1Q24) ,Y4(1024) 

1330  MAT  X4-ZER 
1340  MAT  Y4-ZER 
1350  MAT  X4— X<  J : Ml > 

1360  MAT  Y4-Y < J  s  Ml ) 

1370  Cl-DOT (X4,X4) 

1380  C2-DOT <Y4,Y4) 
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1390  C3=D0T(X4,Y4> 

1400  C4»CiA2 
1410  C5»C2A2 
1420  C6-SUM<X4> 

1430  C7-SUM<Y4> 

1440  N2“N0*C3-C6*C7 

1450  D2«SQR<<N0*C1-C6A2)*<N0*C2-C7A2>> 
1460  R-N2/D2  !  CORRELATION  FACTOR 
1470  R2-RA2 

1480  PRINT  “  CORRELATION  FACTOR  - 
1490  PRINT  “RA23S“}R2 
1500  RETURN 


**“«?«•  ■■i.vf*-*  j*. 


10  •  REDUC2 
20  I 

30  >  BERT  HANSEN 

40  !  THIS  PROGRAM  WILL  RETREIVE  THE  UNIVERSAL  CURVE  DA 
TA  FROM  THE  DISC. 

50  •  IT  WILL  THEN  PLOT  THE  CURVE  AND  THEN  RETREIVE  AND 
PLOT  THE  CALIBRATION  CURVES. 

60  I 

70  CLEAR 
30  OPTION  BASE  1 

90  COM  Q3*li2] ,0241121 ,T1*12Q] ,D1*[121 ,D2*tl2] ,h 

100  D9*»“»D70iM 

110  L-. 000000488 

120  F*S26 

130  K 1S1 . 5365 

140  PLOTTER  IS  705 

150  CSIZE  3, .4,0 

160  OUTPUT  709  ) "TD" 

170  ENTER  709  }  TiO 
180  SHORT  X<100),Y<1QO > 

190  SHORT  T ( 1 024 ) 1 1 ( 1 024  > 

200  P40«"THETA  BAR" 

210  PS*«"INTENS1TY" 

220  DIM  P34C261 

230  P30»" INTENSITY  VS.  THETA  BAR" 

240  LIMIT  0 ,250 ,  0 , 180 
260  LOCATE  20 >120, 20 ,80 
270  FRAME 

280  SCALE  0,10,0,3 
290  LINETYPE  1 
300  E*1 
310  E3-3 

320  E19-"A,.3D"  8  E2f»".3D,A" 

330  Y«i  t  Y2»10 
340  G08UB  1230  !  LABEL 
350  G03UB  1300  «  LOG  SCALE 
360  Eif-"A,.2D"  8  E2*»".2D,A" 

370  E*10*E 

380  Y-20  8  Y2-100 

390  GOSUB  1300 

400  Ei*-"A,D.D"  8  E24-"D.D,A" 

410  Y»200  8  Y2-1000 

420  E-10OE 

430  GOSUB  1391 

440  DEG 

450  LDIR  -91 

460  FOR  X*0  TO  10 

470  •  YAXIS  X 

480  MOVE  X,-.l 

490  LABEL  USING  "DD"  ;  X 
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500  NEXT  X 
510  LDIR  0 

520  DISP  "SCALE  COMPLETE" 

530  SCALE  0,10,0  ,1 
540  CSIZE  4,. 5,0 
550  FXD  1,2 
560  MOVE  5,-. IS 
570  LABEL  P4* 

580  DEG 
590  LDIR  91 
600  MOVE  -1, .5 
610  LABEL  PS* 

620  LDIR  0 
630  MOVE  7, .85 
640  >  LABEL  P3* 

650  MOVE  7,. 95 

660  »  LABEL  "REDUC2  RESULTS" 

670  MOVE  7,. 75 
680  !  LABEL  Tl* 

690  QS*«"UNICUR" 

700  DISP  "  LABELING  COMPLETE" 

710  DISP  "CHANGE  PENS  AND  ICONT1" 

720  PAUSE 

730  SCALE  0,10, -3,0 

740  HS»1024 

750  ASSIGN*  1  TO  85* 

760  FOR  I-i  TO  99 
770  READ*  1  ;  X(I),Y(I) 

780  PLOT  X  < I ) , LGT ( Y  < I ) ) 

790  NEXT  I 
800  PENUP 

810  ASSIGN*  1  TO  * 

820  DISP  "  DO  YOU  WANT  TO  PLOT  A  CALIBRATION  CURVE?  <Y/ 
N>" 

830  INPUT  NO 

840  IF  N*«"N"  THEN  GOSUB  930 

850  DISP  "  ENTER  THE  FILE  NAME  OF  THE  CALIBRATION  CURVE 
TO  BE  PLOTTED.  (CALI, ETC)" 

860  DISP  *  INCLUDE  THE  DISC  DRIVE  WHERE  CAL  FILE  IS  STO 
RED— E.G.  ~NAME s  D701  * 

870  INPUT  Q3t 

880  I  DISP  "ENTER  THE  NUMBER  OF  DATA  POINTS  IN  THE  FILE 
.  (100, 1024, ETC)" 

890  »  INPUT  M2  !  CAL  FILES  HAVE  100  ELEMENTS 
900  M2-100 

910  COSUB  980  i  CALC  THETA  BAR  AND  PLOT  ! 

920  GOTO  820 
930  LINETYPE  1 
940  CLEAR 
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950  DISP  “PLOTS  COMPLETE  •“ 

960  DISP  “LOAD  REDUC3  TO  COMPUTE  D32  FROM  A  DATA  SET .  ** 
970  END 

980  I  READ  AND  PLOT  CAL 

990  DISP  "ENTER  THE  INDEX  OF  REFRACTION  USED  FOR  THE  CA 
L  RUN.  (AIR-1. 0,  UATER-1.3S>“ 

1000  INPUT  H 

1010  DISP  “ENTER  THE  Dbor  USED  TO  COMPUTE  THETA  BAR .  “ 
1020  INPUT  D 

1030  DISP  "ENTER  THE  DIODE  NUMBERS  FOR  THE  REGION  OF  IN 
TEREST.  (Pt ,P2>  TYPICALLY  (SO ,400)“ 

1040  INPUT  Pi,P2 

10SO  H9-D*PI*.02S/(L*F*M> 

1060  Pl«H9*(6i . 46+P1) 

1070  P2«H9*(6i.46-M>2> 

1080  PRINT  “LOWER  LIMIT  OF  THETA  BAR-“iPl 
1090  PRINT 

1100  PRINT  “UPPER  LIMIT  OF  THETA  BAR»“{P2 
1110  ASSIGNO  1  TO  Q3* 

1120  FOR  1*1  TO  M2 
1130  READO  1  }  T(I),I<I> 

1140  NEXT  1 

1150  MOVE  T ( 1 ) , I ( 1 ) 

1160  FOR  1*1  TO  100 

1170  IF  T( I) (Pi  THEN  1200 

1180  IF  T( I ) >P2  THEN  1200 

1190  IF  I(I)>0  THEN  PLOT  T ( I > ,LGT( I ( I > > 

1200  NEXT  I 
1210  PENUP 
1220  RETURN 
1230  !  LABEL 
1240  Y3-Y/1000 

12S0  MOVE  10.1 >LGT( Y>  «  LORG  2 
1260  LABEL  USING  Elf  j  “-“,Y3 
1270  MOVE  -. 1 ,LGT(Y)  §  LORG  8 
1280  LABEL  USING  E20  ;  Y3,“-* 

1290  RETURN 

1300  FOR  Y-Y  TO  Y2  STEP  E 
1310  E2-LGT ( Y  > 

1320  !  XAXIS  LGT(Y) 

1330  IF  RMD<Y,1Q*E>*0  THEN  1350 
1340  GOSUB  1230  !  LABEL 
1350  NEXT  Y 
1360  RETURN 
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10  I  REDUC3 
40  I 

SO  I  THIS  ROUTINE  CONFUTE  THE  VALUE  OF  D32  FOR  THE  DA 
TA  IN  THE  FILE  "8#" 

60  i  BERT  HANSEN  13  JULY  82 
70  I 

80  OPTION  BASE  1 
90  US- 10 24 

100  COH  8301121 ,82#  I 121,T1#[20]>D1#112J ,D2#t 121  ,M 

110  CLEAR 

120  D9#“" :D701" 

130  SHORT  R(20),X1(100) ,Y1 < 100  > 

140  SHORT  X< 1024) , Y< 1024) ,D3<20) ,T3<20) ,T< 1024) 

ISO  DISP  "ENTER  THE  FILENAHE  OF  THE  DATA  TO  BE  REDUCED 
<DATA1:D701,ETO" 

160  INPUT  Q30 

170  ASSIGN#  2  TO  Q30 

180  FOR  I»1  TO  HS-3  STEP  4 

190  READ#  2  i  X(I),Y(I),X(I+1)>Y<I+1),X<1+2),Y(I+2)#X(1 
♦3) > Y< 1+3) 

200  NEXT  I 

210  !  FOR  1-1  TO  1024 
220  I  X(1)-X<I)/198.S6 
230  (  NEXT  I 

240  I  THETA  AND  INTENSITY  ENTERRED 
2S0  La . 000000488 
260  Kl-L/PI 

270  DISP  "ENTER  THE  NAME  OF  THE  DATA  FILE  OF  THE  REFERS 
NCE  ARRAY . <  UNICUR , CAL 1 , CAL2 , CAL 3 ) " 

280  INPUT  Ql# 

290  H4-100 

300  ASSIGN#  1  TO  Ql# 

310  FOR  1-1  TO  H4-1  STEP  2 

320  READ#  1  }  X1<I)>YI(I),X1(I-H)>Y1<I+1) 

330  NEXT  I 

340  DATA  .S,.3,.2>.1,.09,.Q8,'Q7,.Q6,.OS>.Q4,.03>.Q2,.0 

1,  .008, .005 

3S0  FOR  N«1  TO  IS 

360  READ  R(N)  I  REF  INTENSITIES 

370  NEXT  N 

380  DISP  "ENTER  THE  NUMBER  OF  DIODES  TO  BE  ONHITTED  CLO 
8E  TO  THE  CENTERLINE." 

390  INPUT  J 

400  I  D32  CALCULATION 

410  Jl*2 

420  FOR  D*1  TO  IS 

430  R1*R <D)  I  REF  INTENSITY 

440  01-AMIN<Y)  »  SMALLEST  INTENSITY 

450  IF  Rl<01  THEN  GOTO  12S0 

460  FOR  I-Jl  TO  100 
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470  IF  Yi  <  I )  >R1  THEN  GOTO  610 
490  IF  Yl< I )  >-Yl < I-i )  THEN  GOTO  600 
490  I  FIND  THETA  BAR 
500  H2«Y1<I-1)-Y1(I) 

510  K2-X1 < I >-Xi <1-1 ) 

520  H3-Yl<I-t>-Ri 
530  K3-K2*<H3/H2> 

540  T3<D>-Xi<I-l>«-K3  t  THETA  BAR 

550  PRINT  "I-"jI;Yi(I>;Rl»D 

560  PRINT  "I-"}IiT3<D>»Xl<I-l)+K3 

570  PRINT 

580  Ji-I 

590  GOTO  620 

600  PRINT  a  ZERO  OR  POS  SLOPE" 

610  NEXT  I 

620  i  FIND  THE  INTENSITY  CORRESSPONDING  TO  DE8IRED  INT 

ENSITY 

630  M9-1Q00 

640  L9-1 

650  FOR  I-J  TO  H9 
660  IF  Y< I >  >0  THEN  691 

670  PRINT  "INTENSITY  VECTOR  GOING  NEGATIVE  AT  I-"*! 

680  PRINT  "INCREASING  J-J+S"  •  J-J+S  8  GOTO  630 
690  Z-0 

700  FOR  H-I-S  TO  1+5 
710  Z-Z*Y<H> 

720  NEXT  H 
730  Z».l*2 

740  IF  Z)R1  THEN  GOTO  790 

750  !  Z  IS  LESS  THAN  R1 

760  J-I  !  STARTING  POS  FOR  NEXT  SEARCH 

770  DISP  "  I  FOUND  IT" 

780  GOTO  810 

790  L9— Z  !  LAST  AVERAGE 

800  NEXT  I 

810  f  FIND  THETA 

820  K7-X<J)-X<J-t> 

830  H7-L9-Z 
840  H6-L9-R1 
850  K6-K7*<H6/H7) 

860  T4-X<  J-D+K6  I  THETA 
870  !  CONFUTE  D32 
880  D3<D)-K1*T3(D)/T4  !  D32 
890  PRINT 

900  PRINT  "FOR  INTENSITY- " *R<D> 

910  PRINT  "THETA-  ")T4>"  RADIANS" 

920  PRINT  "THETA  BAR-  ")T3<D> 

930  PRINT  "D32-  * ; D3(D> *1000000 j »  MICRONS" 

940  NEXT  0 
950  D5-0 
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960  FOR  11-1  TO  D 
970  D5"D5+D3( II ) 

980  NEXT  XI 

990  DS»D5/D*1000000 

1000  PRINT  "032  AVERAGE*" jOS j "  MICRONS  " 

1010  PRINT 

1028  0I8P  "00  YOU  WANT  TO  COMPARE  OATA  TO  ANOTHER  REFER 
ENCE  FILE?  <V/N>" 

1030  INPUT  NO 

1040  IF  N*-"N"  THEN  1080 

10SO  DI8P  "ENTER  THE  FILENAME  OF  THE  REF.  DATAFILE  (CAL 
1A,ETC)“ 

1060  INPUT  Qlt 
1070  GOTO  290 

1080  DISP  "00  YOU  WANT  TO  REDUCE  ANOTHER  DATAFILE?  (Y/N 
)" 

1090  INPUT  N* 

4  4  01  It  TUTM  4  4  011 

1110  DISP  "ENTER  THE  FILENAME  OF  THE  DATAFILE  TO  BE  RE 
DUCED.  (DATAl ,ETC) " 

1120  INPUT  Q* 

1130  ASSIGN*  1  TO  Q« 

1140  FOR  I-i  TO  1021  STEP  4 

tlSO  READ*  1  i  X<I>,Y<X),X(I>l>,Y<I+l>,X(I+2),Y(I-*-2>,X< 
1+3) , Y< I+3) 

1160  NEXT  I 

1170  GOTO  400  !  REPEAT  D32  CALCULATIONS 
1180  DISP  "REDUCTIONS  COMPLETED” 

1190  DISP  "WANT  TO  DO  ANOTHER  RUN?" 

1200  INPUT  N* 

1210  IF  N*-"N"  THEN  1230 

1220  GOTO  10 

1230  CHAIN  "REDUC4” 

1240  END 

12S0  I  REF  INTENSITY  IS  LESS  THEN  SMALLEST  INTENSITY 
1260  I  DON'T  SEARCH  JUST  COMPUTE  D32  AVERAGE  TO  THIS  PO 
INT 

1270  D5-0 

1280  FOR  11*1  TO  D-l 
1290  DSaDS+D3< ID 
1300  NEXT  XI 

1310  D5*D5/(D-1>  !  D32  AVERAGE 
1320  GOTO  1000 
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10  I  REDUC4 
20  ! 

30  |  BERT  HANSEN 

40  !  THIS  PROGRAM  WILL  RETREIVE  THE  UNIVERSAL  CURVE  DA 
TA  FROM  THE  DISC. 

SO  I  XT  WILL  THEN  PLOT  THE  UNIVERSAL  CURVE. 

60  !  THE  USER  MAY  THEN  PLOT  INTENSITY  VERSUS  THETA  BAR 
WITH  D32  ENTERRED  FROM  KEYBOARD. 

70  ! 

80  CLEAR 

90  GCLEAR 

100  OPTION  BASE  1 

110  COM  030(121 ,0241 12 1,T10( 20), 010(121  >0201 12  ]  #  M 
120  D90-"«D701" 

130  L-. 000000468 
140  F-S26 
ISO  K1-1.S36S 

160  DISP  "ENTER  THE  PLOTTER  ADDRESS  (CRT-1,  722SB-70S)" 

170  INPUT  P 
180  PLOTTER  IS  P 

190  IF  P-70S  THEN  LIMIT  0, 250,0,180 
200  CSIZE  3, .4,0 
210  OUTPUT  709  j"TD" 

220  ENTER  709  j  T10 
230  SHORT  X(10Q) ,Y(10Q) 

240  8HORT  T(1024) ,I< 1024) 

2S0  P 40- "THETA  BAR" 

260  PSO-" INTENSITY" 

270  DIN  P30I261 

260  P30-" INTENSITY  VS.  THETA  BAR" 

99o  champ 

300  LOCATE  20,120,20,80 
310  FRAME 

320  SCALE  0,10,0,3 
330  LINETYPE  1 
340  E-l 
3S0  E3-3 

360  EiO-"A, . 3D"  •  E24-”.3D,A" 

370  Y-l  6  Y2-10 
380  GOSUB  1260  I  LABEL 
390  GOSUB  1330  I  LOG  SCALE 
400  E10-"A,.2D"  •  E2t-".2D,A" 

410  E-10«E 

420  Y-20  •  Y 2-100 

430  GOSUB  1330 

440  EiO-"A,D.D"  •  E24-"D.D,A" 


470  GOSUB  1330 

480  DEG 

490  LDIR  “91 

S00  FOR  X-0  TO  10 

510  YAXXS  X 

520  HOVE  X,-.i 

S30  LABEL  USING  "DD“  }  X 

540  NEXT  X 

5S0  LDIR  I 

560  DISP  "SCALE  COHPLETE" 

570  SCALE  0,10,0,1 
S80  CSIZE  4 , .5,0 
590  FXD  1,2 
600  HOVE  5, “.IS 
610  LABEL  P4* 

620  DEG 
630  LDIR  91 
640  HOVE  -I , . S 
6S0  LABEL  PS* 

660  LDIR  0 

670  HOVE  7,. 85 

680  LABEL  P30 

690  HOVE  7, .95 

700  LABEL  "REDUC4  RESULTS" 

710  HOVE  7,. 75 
720  LABEL  Tl* 

730  QSO*"UNICUR" 

740  DISP  "  LABELING  COHPLETE" 

750  DISP  "CHANGE  PENS  AND  [CONTI" 

760  PAUSE 

770  SCALE  0,10, -3,0 

780  H5-1024 

790  ASSIGN*  1  TO  Q5» 

800  FOR  1-1  TO  99 
810  READ*  1  ;  X<I),Y<I> 

820  PLOT  X< I ) ,LGT< Y ( I ) ) 

830  NEXT  I 
840  PENUP 

850  AS8IGN*  1  TO  0 

860  DISP  "  DO  YOU  WANT  TO  PLOT  A  DATA  FILE  CURVE?  <Y/ 
N>  " 

870  INPUT  N* 

880  IF  N*-"N"  THEN  GOSUB  970 

890  DISP  "  ENTER  THE  FILE  NAME  OF  THE  DATA  FILE  TO  BE 
PLOTTED.  <NAHE;D701) " 

900  DISP  "  INCLUDE  THE  DISC  DRIVE  WHERE  CAL  FILE  IS  STO 
RED— E.G.  -NAHE:D7Q1" 

910  INPUT  Q3* 

920  H2-1024 

930  GOSUB  1010  !  CALC  THETA  BAR  AND  PLOT  t 


941  IF  P*705  THEN  970 

9S0  GRAPH 

960  COPY 

970  LINETYPE  1 

980  CLEAR 

990  DISP  “PLOTS  COMPLETE  1“ 

1000  END 

1010  I  READ  AND  PLOT  CAL 

1020  DISP  "ENTER  THE  INDEX  OF  REFRACTION  USED  FOR  THE  C 
AL  RUN.  < AIR-1. 0,  UATER-1.35)" 

1030  INPUT  M 

1040  DISP  “ENTER  THE  D32  TO  BE  USED  IN  COMPUTING  THETA 
BAR.  " 

10S0  INPUT  D 

1060  DISP  “ENTER  THE  DIODE  NUMBERS  FOR  THE  REGION  OF  IN 
TEREST.  <P1,P2>  TYPICALLY  (50,400)“ 

1070  INPUT  P1,P2 

1080  H9»D*PI*.025/(L*F*N> 

1090  P3*H90 < 61 . 46+P 1 ) 

1100  P4-H9*<61.46+P2> 

1110  K6«D*PI/L 

1120  PRINT  “LOWER  LIMIT  OF  THETA  BAR»“;P3 
1130  PRINT 

1140  PRINT  “UPPER  LIMIT  OF  THETA  BAR-“;P4 
1150  ASSIGN*  1  TO  Q3t 
1160  FOR  1-1  TO  M2  STEP  4 

1170  READ*  1  i  T(I),I(I),T(I+l),I(I+l),T(I+2),I(I+2),T( 
1+3), I (1+3) 

1180  NEXT  I 

119*  MOVE  T<1) ,1(1) 

1200  FOR  I-Pl  TO  P2 
1210  T(I)-K6*T(I> 

122*  IF  I(I)>0  THEN  PLOT  T< I ) ,LGT< I < I ) ) 

1230  NEXT  I 
1240  PENUP 
1250  RETURN 
1260  !  LABEL 
1270  Y3-Y/1000 

1280  HOVE  10.1, LGT ( Y )  8  LORG  2 
1290  LABEL  USING  Elf  *  "-",Y3 
1300  HOVE  - . 1 ,LGT < Y>  •  LORG  3 
1310  LABEL  U8ING  E20  }  Y3,“-“ 

1320  RETURN 

1330  FOR  Y-Y  TO  Y2  STEP  E 
1340  E2-LGT(Y> 

1350  XAXIS  LCT(Y) 

1360  IF  RMD(Y,10*E)*0  THEN  1380 
1370  GOSUB  1260  !  LABEL 
138*  NEXT  Y 
1390  RETURN 
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